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A  Prototype  Windflow  Modeling  System  for 
Tactical  Weather  Support  Operations 


1.  INTRODUCTION 

This  report  is  the  culmination  of  a  three -year  research  effort  at  AFGL  to 
develop  and  evaluate  a  two-dimensional  <x-y  plane)  surface-layer  windflow  model 
for  complex  terrain.  The  report  begins  by  outlining  the  operational  need  that  exists 
today  for  accurate  specification  of  the  low-level  windfield  over  irregular  topography. 
The  model  physics  and  computer  architecture  are  also  described,  as  well  as  the 
peripheral  programs  used  in  the  system.  A  user's  guide  is  then  presented,  which 
is  geared  toward  model  usage  by  operational  forecasters.  We  conclude  with  the 
results  of  an  actual  field  test  made  with  the  model  at  the  Air  Weather  Service  (AWS) 
detachment  at  Ft.  Devens,  Massachusetts.  The  first  section  contains  an  overview 
of  both  the  military  and  nonmilitary  uses  for  the  windflow  modeling  system,  a 
summary  of  previous  research  efforts  with  the  model,  and  instructions  on  the  best 
utilization  of  the  different  sections  of  this  report  that  are  tailored  to  a  variety  of 
audiences. 

1.1  Operational  Need  for  Low-level  Wind  Information 

The  problem  of  specifying  winds  over  irregular  terrain  embodies  many  physical 
considerations  that  impose  difficulties  on  our  ability  to  describe  the  flow  accurately. 

(Received  for  Publication  29  April  1987) 


1 


The  presence  of  both  natural  and  man-made  obstacles  over  the  land  surface,  the 
existence  of  nonhomogeneous  vegetation  and  soils,  bodies  of  water  such  as  lakes 
and  oceans,  and  the  structure  of  the  terrain  itself,  present  a  myriad  of  interacting 
and  nonlinear  physical  processes  that  are  difficult  to  specify  and  sort  out.  The 
advent  of  increasingly  sophisticated  numerical  models  has  helped  us  in  our  under¬ 
standing  of  which  physical  processes  are  important  in  the  formation  of  slope  winds, 
mountain-valley  circulations,  sea/land  breezes  over  complex  coastal  regions,  and 
the  roles  of  these  processes  in  convective  cloud  formation  and  creation  of  meso- 
scale  circulations  and  boundaries.  These  models  range  in  complexity  from  three- 

dimensional  primitive-equation  (PE)  types  such  as  the  Penn  State/NCAR  model 

1  2 
(Anthes  and  Warner  )  and  the  Colorado  State  models  (Mahrer  and  Pielke;  and 

3  4 

Tripoli  and  Cotton  ),  to  simpler  single-layer  models  such  as  those  of  Danard, 

5  6 

Mass  and  Dempsey,  and  Alpert  et  al. 

The  above-mentioned  models  have  been  exclusively  applied  in  research  studies. 

Operationally -available  numerical  models  have  horizontal  and  vertical  resolutions 

that  are  too  coarse  to  resolve  terrain-induced  mesoscale  systems  such  as  mountain- 

7  . 

valley  circulations  and  slope  flows  (Pielke  ).  Obviously,  the  operational  require¬ 
ments  for  specifying  wind  conditions  over  a  small  region  such  as  an  individual 
mountain/ valley  system  make  the  use  of  models  that  run  on  desktop  computers  an 
attractive  option  for  those  tasked  with  mesoscale  weather  forecasting  duties.  Such 
simple  numerical  models  contain  less  physics  than  the  larger,  more  complex  PE 
models.  This  fact  requires  the  forecaster  to  be  more  aware  of  the  model  limita¬ 
tions  and  makes  the  interpretation  of  the  output  a  more  difficult  task. 

1.  Anthes,  R.  A.  ,  and  Warner,  T.T.  ( 1978)  The  development  of  mesoscale  models 

suitable  for  air  pollution  and  other  mesometeorological  studies, 

Mon.  Wea.  Rev.  106:1045-107  8. 

2.  Mahrer,  Y.  ,  and  Pielke,  R.  A.  (1977)  The  effects  of  topography  on  sea  and  land 

breezes  in  a  two-dimensional  numerical  model,  Mon.  Wea.  Rev. 

105:1151-1162.  ~ 

3.  Tripoli,  G..).,  and  Cotton,  W.  R.  (1982)  The  Colorado  State  University  three- 

dimensional  cloud/mesoscale  1982.  Part  I:  General  theoretical  framework 
and  sensitivity  experiments,  J.  Rech.  Atmos.  16:185 -2 19. 

4.  Danard,  M.  (  1977)  A  simple  model  for  mesoscale  effects  of  topography  on 

surface  winds,  Mon.  Wea.  Rev.  105:572-581. 

5.  Mass,  C.  F.  ,  and  Dempsey,  D.  P.  (1985)  A  one -level,  mesoscale  model  for 

diagnosing  surface  winds  in  mountainous  and  coastal  regions,  Mon.  Wea.  Rev. 
113:1211-1227. 

6.  Alpert,  P.  ,  Eppel,  A.,  and  Getenio,  B.  ( 1985)  Surface  wind  prediction  over 

complex  terrain  -  Application  of  a  one-level  terrain  following  model  to 
Israel,  Preprints  Seventh  Conf.  Numerical  Weather  Prediction,  Am. 

Meteorol.  Soc.  ,  Boston,  Mass.,  pp.  369-373. 

7.  Pielke,  R.A.  (1982)  The  role  of  mesoscale  numerical  models  in  very-short- 

range  forecasting,  in  Nowcasting,  K.  A.  Browning,  Ed.,  Academic  Press  Inc., 
New  York. 
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The  AFGL  windflow  modeling  system  described  in  this  report  was  originally 
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designed  by  Ball  and  Johnson  for  use  by  the  Army  Atmospheric  Sciences  Labora¬ 
tory  (ASL)  at  White  Sands,  New  Mexico,  to  specify  terrain-induced  wind  effects 
over  small  areas  of  complex  terrain  (5  to  20  km  on  a  side),  using  high  horizontal 
resolutions  typically  between  100  and  200  m.  The  numerous  mission  support  uses 
for  such  a  model  are  described  below: 

(1)  Low-level  aviation— including  close-air  support  missions,  nap-of- 
the -earth  helicopter  flying,  and  low-level  radar -avoidance  flying 
by  strategic  forces.  Other  missions  are  low-level  cargo  delivery 

by  transport  aircraft,  and  approaches  into  and  out  of  paradrop  zones. 
Diagnosis  of  terrain -induced  turbulence  or  wind  effects  in  or  around 
airfield  locations  is  another  potential  use. 

(2)  Paradrop  operations —a  critical  area  that  is  very  sensitive  to  strong 
or  gusty  surface  winds  and  low-level  wind  shear.  The  setup  of  marker 
panels  over  the  drop  zone  for  the  aircraft  and  jumpers  is  sensitive  to 
wind  direction  and  speed  throughout  the  lowest  300  meters. 

(3)  Windflow  climatology —useful  for  long-  or  short-range  mission  planning. 
This  information  is  critical  over  data-denied  or  data-sparse  areas. 

All  aspects  of  mission  planning  such  as  target  determinations,  locations, 
favorable/unfavorable  times,  and  go  or  no-go  decisions  are  sensitive 
to  this  information. 

(4)  Target -acquisition  winds —important  for  determining  the  types  of 
weapons  to  use  ("smart"  and  conventional).  The  model  is  especially 
useful  for  defining  conditions  on  the  enemy  side  of  the  forward  edge 
of  battle  area  (FEBA)  and  forward-line-of-own  troops  (FLOT). 

(5)  Nuclear,  Biological,  Chemical  (NBC)  Operations— this  may  be  the 
most  wind-sensitive  operation  of  all  since  the  transport  and  diffusion 
of  contaminants  can  be  affected  by  all  aspects  of  the  terrain  and  land 
surface.  It  is  difficult  to  predict  windflow  for  light  wind  situations 

in  which  the  small-scale  surface  features  such  as  vegetation,  soil 
moisture,  and  slope  aspect  are  important  in  producing  certain  flow 
patterns. 

A  considerable  variety  of  nonmilitary  uses  for  this  type  of  model  also  exists. 
Several  of  these  areas  overlap  or  are  complementary  to  the  military  uses  described 
above,  and  thus  need  to  be  mentioned  here: 


8.  Ball,  .1.  A.  ,  and  Johnson,  S.  A.  (1978)  Physically  Based  High  Resolution  Surface 
Wind  and  Temperature  Analysis  for  FTPAMS,  AsL-CR-75^-0043 -1,  U.  Army 
Atmospheric'  Sciences  Laboratory,  White  Sands  Missile  Range,  New  Mexico, 
ADA  055861. 
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(1)  Transport  and  diffusion— obvious  relationship  to  the  NBC  applications, 

especially  for  specifying  hazard  areas  resulting  from  toxic  chemical 

releases.  Another  area  that  has  a  military  complement  is  fire 

weather,  especially  with  regard  to  transport  and  diffusion  of  smoke, 

9 

and  spread  of  the  fire  (see  Svejkovsky  for  description  of  long-range 
smoke  transport  as  detected  from  satellite). 

(2)  Agricultural  meteorology— several  applications  in  this  area  are  frost 
prediction  (inferred  from  examining  cold-air  drainage  and  pooling 
patterns),  canopy  effects  on  the  windflow  (see  Cionco^’  **  for  descrip¬ 
tions  of  a  canopy  parameterization  for  the  ASL  version  of  the  windflow 
model),  soil  erosion  due  to  surface  wind  effects,  and  aerial  spraying 
operations. 

(3)  Wind  energy  siting— specification  of  terrain-induced  effects  on  the 

wind  are  important  for  locating  wind  energy  sites  since  the  windpower 

1 2 

is  proportional  to  the  cube  of  the  wind  speed  (Panofsky  and  Dutton  ). 

Both  numerical  and  empirical  models  have  been  used  to  examine  this 
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problem  (Kirchoff  and  Kaminsky;  Barnard  et  al  ). 

(4)  Civil  aviation  — low  level  aviation  is  sensitive  to  terrain  effects  on  the 
wind,  especially  for  gliders,  light  airplanes  and  helicopters. 

(5)  Numerical  Weather  Prediction  (NWP) — models  such  as  ours  can  be 
used  in  comparison  studies  with  other,  more  complex  models  to 
explore  the  dominant  physical  processes  present  in  terrain-induced  flows. 


9.  Svejkovsky,  .1.  (1985)  Santa  Ana  airflow  observed  from  wildfire  smoke  patterns 
in  satellite  imagery,  Mon.  Wea.  Rev.  113:902-906. 

10.  C'ionco,  R.  M.  (  1983)  On  the  coupling  of  canopy  flow  to  ambient  flow  for  a 
variety  of  vegetation  types  and  densities.  Boundary  -  Layer  Meteorol. 
26:325-335.  '  ' 

1!.  Cionco,  R.  M.  (1985)  On  modeling  canopy  flow  coupled  to  the  surface  boundary 
layer,  Proc.  17th  C'onf.  on  Agriculture  and  Forest  Meteorology  and  7th  Conf. 
on  Biometeorologv  and  Aerobiology,  Am.  Meteorol.  Soc.  ,  Boston,  Mass.  , 
pp.  116-119. 

12.  Panofskv,  11.  A.  ,  and  Dutton,  .LA.  (1984)  Atmospheric  Turbulence;  Models 

and  Methods  for  Engineering  Applications,  Wiley,  New  York,  297  p. 

13.  Kirchoff,  R.H.,  and  Kaminsky,  F.  C.  (1983)  Empirical  Modeling  of  Wind 

Speed  profiles  in  Complex  Terrain,  DOE  -ET  /  1037  4-  82  /  1  (DE  839  1  1 6  13 ), 
Pacific  Northwest  Laboratory  (Contract  1)1  -AC06-76RLO  1830), 

Richland,  Wash.  ,  99352. 

14.  Barnard,  .I.C.,  Weglev,  H.  L.  ,  and  Hiester,  T.  R.  (1985)  Improving  the 

Performance  of  Mass -Cons istent  Numerical  Models  Using  Optimization 
Techniques,  PNl. -’566  (Contract  DE -AC06-?6f?LO  1830),  Pacific  Northwest 
Laboratory,  Richland,  Wash.  ,  99352, 


4 


A  windflow  model  analysis  can  be  used  to  initialize  a  high-resolution 

PE  model  for  areas  of  complex  terrain  where  an  objective  analysis 

is  unsuitable.  Prediction  of  other  weather  phenomena  such  as 

radiation  fog  could  be  accomplished  by  examining  cold -air  drainage 

patterns  to  infer  or  actually  calculate  temperature  and  moisture 

15 

distributions  (Tabony  ). 

1.2  Previous  Work  With  the  Model 

The  research  and  development  involving  this  model  over  the  last  eight  years 
has  mainly  centered  on  parallel  efforts  at  ASL  and  AFGL. 

At  ASL,  the  model  has  been  employed  in  conjunction  with  a  pseudo -particle 

diffusion  model  to  study  the  dispersion  of  smoke  in  battlefield  environments 

16  17  18 

(Cionco;  Ohmstede;  and  Ohmstede  and  Stenmark  ).  Further  development  of 

19 

the  windflow  model  was  done  by  Amlicke  and  Coleman,  who  examined  the  way 

the  model  determines  the  exponent  for  the  power -law  wind  profile,  and  created  a 

capability  for  the  model  to  be  initialized  with  the  modified  output  from  a  larger- 

scale  numerical  model.  Cionco^’  **  described  the  coupling  of  a  vegetative  canopy 

formulation  to  the  model,  and  documented  the  windflow  model' s  place  in  ASL's 

20 

model  hierarchy  (Cionco  ).  The  ASL  model  hierarchy  is  a  system  of  numerical 

models  that  covers  the  mesoscale  from  horizontal  scales  of  200  km  down  to  5  km. 
2 1 

Veazey  and  Tabor  used  the  windflow  model  to  examine  requirements  for  meteor¬ 
ological  sensor  density  on  the  battlefield. 

15.  Tabony,  R.C.  (1985)  Relations  between  minimum  temperature  and  topography 

in  Great  Britain,  J.  Climatology  J3:503-520. 

16.  Cionco,  R.M.  (1982)  A  meteorological  approach  to  chemical  defense  over 

complex  terrain  with  vegetation.  Workshop  on  the  Parameterization  of 
Mixed-Layer  Diffusion,  Las  Cruces,  N.  Mex.  ,  pp.  323-328. 

17.  Ohmstede,  W.  D.  (1982)  The  parameterization  of  battlefield  dispersion-new 

frontiers.  Workshop  on  the  Parameterization  of  Mixed-Layer  Diffusion, 

Las  Cruces,  N.  Mex.,  pp.  279-287. 

18.  Ohmstede,  W.  D. ,  and  Stenmark,  E.  B.  (1982)  A  model  for  characterizing 

transport  and  diffusion  of  air  pollution  in  the  battlefield  environment, 
Workshop  on  the  Parameterization  of  Mixed-Layer  Diffusion,  Las  Cruces, 

N.  Mex.,  pp.  416-423. 

19.  Amlicke,  B.B.,  and  Coleman,  I.W.  (1984)  High  Resolution  Wind  (HRW)  Model, 

MRC/ WDR -R -089,  Mission  Research  Corporation,  Alexandria,  VA. 

20.  Cionco,  R.M.  ( 1985)  Modeling  airflow  over  variable  terrain,  Proc.  of  the 

HAZMAT  '85  West  Conference,  Long  Beach,  Calif. 

21.  Veazey,  D.R.,  and  Tabor,  P.  A.  ( 1985)  Meteorological  sensor  density  on  the 

battlefield,  Workshop  on  Geographic  Information  Systems  in  Government, 
Bruce  Opitz  (Ed.  ),  Deepak  Publishing,  Hampton,  VA,  pp.  195-208. 


The  model  was  acquired  by  AFGL  in  1984  for  incorporation  into  a  toxic  chemi- 
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cal  prediction  system  for  complex  terrain  (Weber  ).  It  has  been  used  to  study  the 

sensitivity  of  surface -layer  flow  to  the  hori7ontal  distribution  of  vegetative  cover 

under  various  stability  conditions  (Lanicci  ).  A  real-data  study  was  performed 

on  the  model  using  meteorological  tower  data  from  Vandenberg  AFB,  California 
24 

(Lanicci  and  Weber  ).  Further  model  developments  at  AFGL  included  the  creation 
of  different  meteorological  data  input  options  (including  a  cursory  objective  analysis 
for  multiple  observations),  and  the  adaptation  of  the  code  to  small  computers  such 
as  the  Zenith -100.  The  model  has  been  applied  to  various  types  of  terrain,  includ¬ 
ing  gentle  topography  with  complex  vegetation  cover  (Ft.  Polk,  Louisiana),  rough 
terrain  with  sparse  vegetation  (Vandenberg),  and  rolling  terrain  with  fairly  uni¬ 
form  forest  cover  (Ft.  Devens,  Massachusetts). 

1.3  How  to  Use  This  Report 

We  have  designed  this  report  to  be  as  comprehensive  as  possible  in  the  docu¬ 
mentation  of  the  modeling  system.  Since  it  is  nearly  impossible  to  present  a  text 
that  is  universally  applicable  to  operational  and  research  people  alike,  we  have 
designated  separate  sections  aimed  at  different  readers. 

The  AFGL  windflow  modeling  system  is  available  for  varied  groups  to  use  for 
the  many  applications  described  in  Section  1.  1.  Researchers  will  find  Sections  1, 

2,  3  and  possibly  Section  5  to  be  most  useful  to  them.  Computer  scientists  who  are 
primarily  programmers  will  most  likely  have  use  for  the  flowcharts  and  subroutine 
descriptions  in  Section  3,  and  find  Section  4  useful  if  they  do  not  have  a  meteorolo¬ 
gical  background  but  wish  to  understand  how  the  model  works.  Operational  users 
such  as  weather  forecasters  will  find  Section  1.  1  to  be  helpful,  as  well  as  Sections 
4  and  5,  and  the  descriptions  of  the  Ft.  Devens  case  studies  in  Appendix  A.  Any¬ 
one  wanting  to  develop  this  system  further,  or  use  the  system  for  their  own  opera¬ 
tional  purposes  will  want  to  read  this  entire  report. 

22.  Weber,  H.  (1986)  private  communication. 

23.  Lanicci,  J.  M.  (1985)  Sensitivity  Tests  of  a  Surface -Layer  Windflow  Model  to 

Effects  of  Stability  and  Vegetation,  AFGL-TR -85-0265,  ADA  169136. 

24.  Lanicci,  J.  M.  ,  and  Weber,  H.  (1986)  Validation  of  a  Surface-Layer  WindFIow 

Model  Using  Climatology  and  Meteorological  Tower  bata  1'roni 


The  guidelines  on  the  preceding  page  are  just  suggestions.  We  hope  they  are 
useful  to  those  readers  who  may  be  unfamiliar  with  models  of  this  type  and  need 
some  direction  in  finding  out  more  about  our  system  in  the  least  confusing  manner. 


2.  THE  MODEL  AND  ITS  PHYSICS 

The  basic  theory  and  equations  for  the  windflow  model  are  contained  in  Ball 
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and  Johnson,  Amlicke  and  Coleman,  and  Lanicci.  We  present  here  a  summary 
of  the  theoretical  framework  surrounding  the  model  equations  and  an  outline  of  the 
model's  iterative  procedure  to  arrive  at  its  windflow  analysis. 

2.1  Variational  Analysis  Theory 

The  theoretical  basis  for  the  model  equation  is  Gauss'  Principle  of  Least 
Constraint.  The  essence  of  the  theory  is  that  the  motion  that  occurs  in  nature 
takes  place  in  such  a  way  as  to  minimize  constraint  forces  arising  from  kinematic 
conditions.  In  order  to  better  understand  Gauss'  principle,  we  proceed  through 
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the  following  analyses  taken  from  Lanczos,  and  Ball  and  Johnson. 

2.  1.  1  FOR  A  SYSTEM  OF  POINT  PARTICLES 

We  begin  with  D'Alembert's  Principle  of  Virtual  Work,  which  states  that  any 
system  of  forces  is  in  equilibrium  if  we  add  the  inertial  (created  by  the  motion) 
forces  to  the  impressed  forces: 

N 

^  (Fk  "  mk^k)  '  6  Rk  =  0  *  (1) 

k  =  1 

where  the  summation  is  taken  over  the  system  of  N  particles,  Ffe  •  fiR^  is  the 
virtual  work  done  on  a  system  by  the  impressed  forces,  and  -m^A^  is  the 
inertial  force. 

Gauss  applied  D'Alembert's  principle  at  time  t  +  r  where  T  is  an  arbitrarily 
small  time  interval,  thus  choosing  5  R  (t  +  t)  to  be  embodied  as  a  variation  in  the 
acceleration  A  : 

5R(t+T)  *  1/2  5A  (t)  •  t2  .  (2) 


25.  Lanczos,  C.  ( 1970)  The  Variational  Principles  of  Mechanics  (4th  Ed.  ) 
U.  of  Toronto  Press,  Toronto,  pp.  106-110. 


which  is  derived  from  Taylor's  expansion  of  the  particle  path  deviation  5  R^.  We 
rewrite  Eq.  (1)  as: 

N 


S  (Fk  -  m^A^)  •  6^  =  0 

k_  i  -  -  K 


(3) 


Since  is  given  and  cannot  be  varied,  we  can  rewrite  Eq.  (3)  as  shown  below: 


y*  —  -*  (F.  -  A^l 

(F,  -  m,  A,  )  •  5  — - - -  =  0  , 

k  =  1  K 


k  k 


m 


(4) 


then  applying  the  derivative  rule  for  the  quantity  in  Eq.  (4)  yields: 


N 

6  >: 


—  (F.  -  m.  A.  )2  =  0 
k  _  j  £.mk  k  k  k 


(5) 


Gauss  redefined  Eq.  (5)  as  follows: 


N 

z  =  D 

k  -  1 


sr-  <l?k  '  mk\'2 


(6) 


where  the  quantity  Z  is  defined  as  the  "constraint"  of  the  motion.  The  principle 
in  Eq.  (6)  says  that  motion  occurring  in  nature  seeks  to  minimize  Z.  If  there  are 
no  constraints,  Z  =  0,  and  Eq.  < 6 >  is  reduced  to  Newton's  Law  of  Motion.  This 
principle  is  exactly  analogous  to  the  Principle  of  Least  Squares  or  Least  Action 
that  forms  the  basis  for  many  variational  analysis  schemes  that  seek  to  mini¬ 


mize  the  difference  between  a  theoretical  value  of  a  function  and  its  observed  value. 
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See  the  discussions  bv  Sasaki,  ’  and  Haltiner  and  Williams,  for  examples  of 
the  application  of  (lie  Least  Squares  principle  to  numerical  analysis. 


26.  Sasaki,  V.  (1058)  An  objective  analysis  based  on  the  variational  method, 

.1.  Meteorol.  Soc.  Japan  36:77  -88. 

27.  Sasaki,  V.  (1070)  Some  basic  formalisms  in  numerical  variational  analysis, 

Mon,  Wca,  Rev.  ^8:87.6-883. 

28.  Haltiner,  G.J.,  and  Williams,  K.T.  (1080)  Numerical  Prediction  and  Dynamic 

Meteorology  (2nd  Ed.  ),  Wiley  and  Sons,  New  York,  447  p. 


2.  1.  2  FLUID  MECHANICS  ANALOGY 


Within  the  realm  of  fluid  mechanics,  there  have  been  several  applications  of 
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a  least  action  principle  approach  to  describing  fluid  motion.  Dutton  presented 
a  derivation  of  the  meteorological  equation  of  motion  using  a  least  action  approach. 
The  constraints  imposed  upon  the  system  were  isentropic  motion  (an  energy  con¬ 
straint)  and  inviscid  conditions.  We  present  the  results  of  his  derivation  below 
[his  Eq.  (16)  in  Chapter  11,  Section  5J* 

AA  *  -  J  J  p5x  •  (  +  v  •  Vv  +  i-7p  +  7$)  dV  dt 

"  If  PT  5s  dV  dt  ,  (7) 

where  AA  is  the  variation  of  the  action  A  of  the  system,  p  is  the  density,  v  is  the 

wind  velocity  vector,  p  is  pressure,  ♦is  the  geopotential  (potential  energy  of  a 

unit  mass  in  the  earth's  gravity  field),  T  is  the  temperature,  and  6s  is  the 

variation  of  specific  entropy  s.  Equation  (7)  contains  two  integrals  over  the  fluid 

volume  and  time,  and  Dutton  states  that  for  isentropic  motion,  6s  *  0,  thus  leaving  only 

the  first  integral.  For  the  action  A  to  be  a  minimum.  A  A  must  vanish  for  arbitrary 

variation  6x  vanishing  outside  the  time  interval  t  ,  t _ ,  and  on  the  fluid  boundaries. 

o  1 

This  condition  necessitates  that  the  quantity  in  parentheses  in  Eq.  (7)  vanish,  which 
is  the  equation  of  motion.  Dutton's  derivation  demonstrates,  from  an  energetics 
approach,  that  the  atmospheric  kinetic  energy  is  only  as  large  as  is  needed  to 
respond  to  external  forces  acting  upon  it. 

g 

Ball  and  Johnson  took  another  approach  in  describing  the  fluid  motion  by 
applying  Gauss'  Principle  of  Least  Constraint  to  incompressible,  inviscid  flows: 

6  f  (A  +  g)“  dV  -  6  f  t  A  •  ds  *  9  .  (8) 

v  s 

where  A  represents  the  fluid  acceleration,  g  is  gravity  (an  external  force),  and  T 
is  he  scalar  stress  in  the  outward  normal  direction  to  the  surface  s.  These  two 
integrals  are  taken  over  the  fluid  volume  and  its  boundary  surface*,  respectively. 
Notice  the  similarity  between  the  volume  integral  in  Eq.  (8)  and  the  formulation  for 
point  particles  shown  in  Eq.  (5).  Bv  applying  the  incompressibility  assumption 
through  introduction  of  Lagrange  multipliers  X  and  applying  some  manipulations  to 
Eq.  (81,  Ball  and  Johnson  produce  a  form  of  the  Gaussian  principle  shown  below 
[  see  their  Eqs,  (11.2)  and  following!: 


(9) 


f( A  +  g  +  -  7X)  ■  6A  dV  -  f  i  (r  +•  X)  5  A  •  ds  =  0 

v  P  a  P 


We  can  expand  Eq.  (9)  and  set  X  =  p  (the  pressure)  to  produce  another  form: 


/  (■££■  +  v  •  Vv  fg  t-  —  Vp)  •  5A  dV  -  J  —  (t  +  p)  6  A  •  ds  -  0.  (10) 


Notice  from  Eq.  (10)  that  the  term  in  the  parentheses  of  the  volume  integral  repre¬ 
sents  the  equation  of  motion.  Notice  that  the  minimum  constraint  presented  here 
corresponds  to  a  situation  where  the  equation  of  motion  and  the  kinematic  constraint 
are  both  satisfied.  The  pressure  gradient  term  in  Eq.  (10)  results  from  the  con¬ 
dition  of  incompressibility,  thus  making  the  pressure  gradient  a  force  of  constraint 
arising  from  incompressibility. 

To  apply  Gauss'  principle  to  the  surface  layer.  Ball  and  Johnson  assumed  that 
the  fluid  volume  is  a  single  layer,  and  introduced  some  simplifying  assumptions  to 
derive  the  basic  windflow  model  equation.  The  first  oT  these  is  the  introduction  of 
a  Boussinesq  approximation  to  characterize  the  external  force  g  as  a  buovancy 
force  B,  defined  in  Eq.  (11): 


B  -  p’g/pr 


(11) 


where  p'  is  the  density  perturbation  from  its  ambient  value  pQ.  This  assumption 
along  with  the  Archimedean  principle,  allows  the  buoyancy  force  to  be  described 
using  the  potential  temperature  0  instead  of  density: 


b  =  sel  =  zgL 
Po  X 


_  (ft  -ft  ) 

o  s 

8  -7) - 


(12) 


where  ft^  is  the  surface  potential  temperature,  and  ft^  is  a  reference  potential 

temperature  above  the  surface.  Since  we  want  to  applv  these  relations  to  the 

surface  laver,  0  should  be  chosen  so  that  Ift  -  ft  )  accurately  represents  the  surface 
o  os' 

layer  stability.  The  Boussinesq  assumptions  can  be  applied  to  Eq.  (8)  with  the 
following  result: 


(ft  -  ft  )' 

A  +  g  -°  7)— 


dV 


ds  =  0  . 


(12) 
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YVe  have  also  introduced  the  dynamic  pressure  p'  in  the  surface  integral  as  an 
additional  external  force.  The  next  simplification  into  the  system  is  the  introduc¬ 
tion  of  terrain-following  coordinates.  This  allows  some  terms  in  Eq.  (13)  to  be 
discarded  in  the  application  to  the  single-layer  representation.  A  no-slip  condi¬ 
tion  at  the  surface  [V(Z  =  0 )  =  0]  means  that  the  only  contribution  to  the  surface 
integral  in  Eq.  (13)  is  from  the  upper  boundary,  which  is  the  boundary  between  the 
surface  layer  and  atmospheric  layers  above.  Since  the  model  only  considers  a 
single  layer,  the  surface  integral  and  any  surface  normal  components  of  the  volume 
integral  are  discarded  as  a  modeling  approximation,  leading  to  the  following 
equation: 


v  o 


dV  =  0  , 


(14) 


where  the  "ll"  denotes  surface  parallel  components.  The  surface -parallel  compo¬ 
nents  of  buovanev  in  Eq.  (14)  result  in  the  introduction  of  horizontal  temperature 

gradients,  which  in  turn  introduce  a  horizontal  pressure  gradient  force,  the  latter 
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force  being  one  of  the  driving  mechanisms  involved  in  slope  flows  (Atkinson  ).  The 
expression  in  Eq.  (14)  can  be  redefined  in  the  same  fashion  as  Eq.  (5)  was,  pro¬ 
ducing  the  basic  model  equation  shown  below: 


R 


t 


f  A 


<e  -  w  n 

o  s 

— 0 - 

o 


dV  , 


(  15) 

- 


where  R  is  the  total  constraint  R  over  the  modeled  volume,  that  is  to  be  minimized 
bv  the  integration  of  Eq.  (15).  One  final  note  about  Eq.  (15)  is  the  assumption  of 

3  V 

steady-state  conditions  (that  is,  -gq-=  0)  in  the  momentum  field.  This  means  that 
only  the  advert  ion  V '  ■  V  V  contributes  to  the  acceleration  term  A.  This  assumption 
means  that  the  windflow  model  is  essentially  an  analysis  tool,  the  implications  of 
winch  are  further  explored  in  Section  2.3. 

2.2  Relaxation  Procedure  and  Physical  Parameterization* 

Having  established  the  basic  model  equation  j  Eq.  (15)],  we  now  describe  the 
means  bv  which  the  minimum  constraint  R  is  reached.  The  version  of  the  model 
described  in  this  report  uses  a  single  surface  observation  of  wind  direction  and 
speed,  temperature,  and  (  loud  cover,  along  with  date  and  time,  to  initialize  the 


30.  Atkinson,  B.  W.  (1981)  Mesoscale  Atmospheric  Circulations,  Academic  Press, 
London,  pp.  252-254. 


model  with  a  uniform  windfield  and  buoyancy.  The  model  produces  its  terrain- 

induced  windflow  analysis  through  a  relaxation  procedure,  incrementally  adjusting 

the  initially  uniform  windfield  at  each  grid  point  every  iteration  step.  The  iterative 

procedure  continues  in  the  model  until  the  first  minimum  of  R  over  the  domain  is 
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achieved.  Numerical  experiments  performed  by  Lanicci  and  Weber  over  artificial 
terrain  showed  that  integration  of  the  model  equation  beyond  the  first  Rt  minimum 
produced  oscillations  in  R{,  eventually  reaching  a  near-absolute  minimum  after 
about  195  steps.  However,  the  resultant  windfield  analysis  is  a  pure  slope-wind 
solution  that  is  independent  of  the  ambient  initial  wind  speeds.  Therefore,  the 
first  or  local  R{  minimum  is  used  as  the  point  at  which  the  model  integration  pro¬ 
cedure  stops,  as  it  produces  a  solution  that  is  meterologically  consistent,  having 
well-integrated  the  effects  of  topography,  stability,  mass  conservation,  and 
momentum  advection.  The  relaxation  procedure  is  described  by  the  flow  chart  in 
Figure  1. 
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Figure  1.  Flowchart  Depicting  Iterative  Procedure  by  Which  Model  Tries  to 
Minimize  Constraint  in  Eq.  (15) 


The  terrain-following  coordinate  system  employed  in  the  model  uses  a  stag¬ 
gered  grid  system.  Terrain  elevations  are  on  one  set  of  grid  points,  with  wind 
calculation  grid  points  staggered  in  between,  so  that  calculated  terrain  slopes  can 
be  used  in  the  model  integration.  The  single-layer  formulation  of  the  model  is 
represented  by  using  the  flux  form  of  the  momentum  advection  in  Eq.  (15).  The 
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model  calculates  momentum  flux  through  boxes  of  thickness  z  ,  whose  corners 
correspond  to  the  terrain  grid  points.  A  schematic  diagram  of  these  flux  box 
elements  on  the  model  domain  is  shown  in  Figure  2.  The  momentum  flux  through 
each  face  of  the  flux  box  is  calculated  and  modified  by  the  effects  of  the  local 
surface-layer  stability  to  determine  the  contribution  of  each  flux  box  element  to 
the  A  term  of  Eq.  (15). 


HORIZONTAL  REFERENCE  PLANE 


Figure  I  Man  rant  Slummy  Grid  Structure'  and  Geometry  of  Model  Domain, 
With  Structure'  of  Individual  Flux  Box.  Notice  that  the  staggered  grid  system 
allows  terrain  slopes  centered  on  the  wind  grid  points  to  be  integrated  into 
model  i  ih  ul.it  ions,  through  the  terrain-following  coordinate  axes  a.  and  a.} 

( a  *  t  e  r  Ball  in  I  Johnson1^) 


I  M<  basis  'fit  i.n  cm  irnriL'  the  surface -layer  stability  is  the  original  calcula- 


•m>ii  oi  ihe-  buoy  mu  v  t> a  n  B,  u-o;ig  the  surface  observation  to  diagnose  a  stability 
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class  based  on  a  Pasquill-type  category  formulation  (Turner  ).  Kunkel  devised 
a  continuous  stability  parameter  SP  for  the  AFGL  diffusion  model  that  we  have 
adapted  for  diagnosing  stability  in  the  windflow  model.  The  stability  parameter  is 
analogous  to  the  Pasquill  category,  but  its  use  avoids  sudden  changes  in  stability 
and  thus  hazard  distance,  with  slight  changes  in  windspeed,  solar  elevation  angle, 
or  cloud  cover.  The  values  of  SP  range  from  0.  5  (unstable -Pasquill  category  A) 
to  6.0  (stable-Pasquill  category  F),  and  their  relationship  is  illustrated  in  Table  1. 

Table  1.  Relation  Between  the  Pasquill  Stability  Categories  (SC)  and  the 
Continuous  Stability  Parameter  (SP).  (After  Kunkel)^ 


SC 

A 

B 

C 

D 

E 

F 

SP 

0.  5-1 

1-2 

2-3 

3-4 

4-5 

5-6 

The  scheme  for  calculating  SP  based  upon  the  surface  observation  is  described 
32 

in  detail  by  Kunkel.  We  summarize  the  method  using  the  flowchart  of  Figure  3. 
The  latitude,  longitude,  date,  and  time  are  used  to  determine  incoming  solar 
radiation  as  a  function  of  solar  elevation  angle,  according  to: 


I  =  S  (t  -  t  )  sin  ( A )  r 
o  1  o 


csc(  A ) 


(16) 


where  I  is  the  incoming  solar  radiation  over  a  unit  horizontal  area  on  the  earth's 

■  2 

surface  during  time  (t,  -t  ),  S  is  the  solar  constant  (1353.3  Wm  ),  A  is  the  solar 

loo 

elevation  angle,  and  r  is  the  atmospheric  transmission  coefficient  (set  equal  to 
0.7  in  the  model).  Also  notice  from  Figure  3  that  the  presence  of  clouds  enters 
into  the  stability  calculation  in  two  wavs.  The  first  way  is  during  the  day,  when  a 
cloud  cover  of  7/8  or  8/8  (overcast  conditions)  causes  the  1  calculation  in  Eq.  (16) 
to  be  modified  bv  the  presence  of  different  types  of  clouds,  using  the  formulations 


ontained  in  Table  17  of  the  Smithsonian  Meteorological  Tables 
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The  second  wav 


31.  Turner,  I).  H.  M'*64)  A  diffusion  model  for  an  urban  area,  .1.  Appl.  Meteorol. 

3:83-"  1. 

32.  Kunkel,  H.  A.  (T'85)  Development  of  Atmospheric  Diffusion  Model  for  Toxic 

Chemical  Releases”  ^ F’GT.  -TR  -85  -0338,  ADA  l6(M3?>. 

3  t.  [.ist,  R.d.,  Ed.  <  l'*83)  Smithsonian  Meteorologi cal  doubles  (\ol.  114) 

Sm  ithsonian  Misc  el  laneous  Collections. 
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Figure  3.  Flowchart  Describing  Procedure  for  Model  Diagnosis  of  Stability 
Parameter  and  Buoyancy  Magnitude,  Using  Input  of  Surface  Observation 


is  at  night,  when  the  cloud  cover  modifies  the  computation  of  surface  heat  flux  H 
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through  the  relation  shown  in  Eq.  (17)  (after  Smith'  ): 


34,  Smith,  F.  B.  (1TT72)  A  scheme  for  estimating  the  vertical  dispersion  of  a  plume 
from  a  source  near  ground  level,  Proc.  of  the  Third  Meeting  of  the  Expert 
Panel  on  Air  Pollution  Modeling,  NATO  Committee  on  the  Challenges  of 
Modern  Society,  Paris,  France. 


where  Hq  =  -40  Wm  (downward  nighttime  heat  flux  in  the  model),  and  c  is  the 
cloud  cover  (varying  from  0  to  8).  During  the  daytime,  the  formula  for  H  becomes: 


H  =  0.4  (I  -  100)  . 


Using  the  I  value  computed  from  Eq.  (16).  An  iterative  procedure  to  calculate 

friction  velocity  u.  and  Monin -Obukhov  length  L  is  then  performed,  using  formulas 

35 

from  Ragland  and  Dennis  based  on  different  magnitudes  of  H.  Once  the  values 

of  ut  and  L  have  been  determined  by  the  appropriate  formulas,  the  stability  param- 

w  36 

eter  SP  is  computed  by  applying  the  nomogram  defined  by  Golder,  as  put  in 
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equation  form  by  Kunkel  : 

SP  =  A  +  B  log  Z  ,  (19) 

where  A  and  B  are  constants  that  are  functions  of  1/L  and  Zq  is  the  surface  rough¬ 
ness  length.  We  should  note  here  that  Zq  is  constant  over  the  model  domain  in 
the  version  described  in  this  report.  If  one  desires  to  use  a  nonhomogeneous  Zq, 
an  appropriate  value  must  be  chosen  for  use  in  Eq.  (19),  or  one  could  apply  Eq.  (19) 
at  everv  grid  point,  using  the  individual  grid  value  of  Zq,  that  would  yield  values 
of  SP  over  domain  as  a  function  of  nonhomogeneous  surface  roughness.  Once  SP 
has  been  determined  by  Eq.  (19),  the  buoyancy  B  can  be  calculated  using  a  form 
of  Eq.  (12)  modified  for  use  of  temperature  instead  of  potential  temperature: 

S„  ■  ?  «0I 


where  AT  is  the  vertical  temperature  change  over  depth  AZ  (100  m),  y  is  the  adia¬ 
batic  lapse  rate  (=  1.0°C  100  m  *),  and  Tg  is  the  surface  temperature  in  K.  The 

vertical  temperature  change  AT  is  a  function  of  SP,  and  they  are  related  to  each 
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other  through  the  values  found  in  U.  S.  NRC  Regulatory  Guide  1.23,  and 


35.  Ragland,  K.  W.  ,  and  Dennis,  B.  L.  ( 1975)  Point  source  atmospheric  diffusion 

model  with  variable  wind  and  diffusivity  profiles,  Atmos.  Environ.  , 

9: 175- 189. 

36.  Golder,  I).  (1972)  Relations  among  stability  parameters  in  the  surface  layer. 

Boundary  Layer  Meteorol.  3:47  -58. 

37.  USNRC’  (1972)  Onsite  Meteorological  Programs  Regulatory  Guide  1.23, 

L.  S.  Nuclear  Regulatory  Commission. 


summarized  by  Table  1  in  Sedefian  and  Bennett.  Figure  4  shows  a  graph  relating 
SP  to  the  actual  lapse  rate  AT/AZ  (solid  curve).  The  model  uses  an  exponential 
equation  fit  to  this  curve,  that  is  described  by  the  dashed  curve  in  Figure  4. 


Figure  4.  Graph  Showing  Stability  Parameter  SP  vs  Vertical  Temperature 
Gradient  (°C  100  m  '  *)  From  NRC  Guide  (solid  line)  and  Model  Equation  Fit 
(dashed  line).  Differences  between  the  two  curves  are  shown  by  stippling. 
Model  overestimates  of  the  magnitude  of  AT/AZ  are  shown  by  light 
stippling,  underestimates  by  dark  stippling 


38.  Sedefian,  L.  ,  and  Bennett,  E.  ( 1980)  A  comparison  of  turbulence  classifica¬ 
tion  schemes,  Atmos.  Environ.  14:741-750. 


The  shaded  areas  show  the  differences  between  the  exponential  formulation  and  the 

NRC  Guide  values  of  AT/AZ  vs  SP.  The  largest  discrepancies  are  for  stability 

parameters  approaching  6.  0,  where  the  model  curve  overestimates  AT/AZ  by  as 

much  as  1°C.  Despite  these  differences,  the  correlation  coefficient  for  the  model 

equation  is  around  0.  99.  More  important  than  these  discrepancies  in  AT/AZ  are 

those  in  buoyancy  magnitude  B  as  a  result  of  the  exponential  curve  fit.  These 

differences  are  highlighted  in  Figure  5,  which  shows  B  as  a  function  of  SP  for 

the  NRC  table  (solid  curve),  and  for  the  exponential  solution  from  Figure  4.  As 

in  Figure  4,  the  largest  discrepancies  appear  for  stable  cases  as  SP  approaches 
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6.  0.  These  differences  are  only  as  large  as  0.  035  ms  ,  and  since  Lanicci 
found  that  the  model  windfields  are  not  as  sensitive  to  the  buoyancy  magnitude  as 
they  are  to  its  sign,  we  believe  that  this  curve  fit  is  adequate  for  diagnosing 
stability  conditions  in  the  model. 

Once  the  buoyancy  has  been  computed,  it  remains  constant  over  the  domain 

throughout  the  entire  integration.  However,  the  local  surface-layer  stability  at 

each  grid  point  is  determined  by  calculating  the  value  of  a  power-law  wind  profile 

exponent.  This  procedure  is  summarized  by  the  upper  flowchart  shown  in  Figure  6. 

The  wind  profile  exponent  n  is  influenced  by  the  surface-layer  stability,  using  the 

12 

equation  below  (after  Panofsky  and  Dutton  ): 

n  =  0/[lnZ/Z)+i//),  (21) 

* m  1  o  ^m 

where  d  is  the  Monin-Obukhov  nondimensional  vertical  wind  shear,  and  \ b 
m  Dm 

is  a  Monin-Obukhov  function  dependent  upon  1/L.  The  values  of  <b  and  li'  are 

39  m  Dn 

determined  through  use  of  the  Businger  et  al  profiles  at  each  grid  point.  The 

value  of  Z  is  determined  by  adding  Zq  to  the  model  calculation  height, 

set  at  10  m  above  ground.  Since  we  are  interested  in  windflow  close  to  the  surface, 

10  m  is  a  good  calculation  height  to  use  since  the  flow  at  this  level  is  often  strongly 

12 

influenced  bv  surface  features  (Panofsky  and  Dutton  ).  The  upper  flowchart  in 

Figure  6  shows  how  the  bulk  Richardson  number  Ri^  determines  the  stability  regime 

used  to  calculate  appropriate  values  of  <t>  ,  1 1/  ,  and  thus  n.  The  Ri,  formula  used 

m  ’m  b 

in  this  model  is  shown  below: 

Rt.  -  BZ/u“  .  (22) 

b 

3'1.  Businger.  .I.A..  Wvngaard,  .1.  C.  ,  Izumi,  Y.  ,  and  Bradley,  E.  F.  (1971)  Flux- 
profile  relationships  in  the  atmospheric  boundary  layer,  J.  Atmos.  Sci. 
28:181-189. 


Figure  5.  Graph  Relating  Stability  Parameter  SP  vs  Buoyancy  Magnitude  for 
NRC  Guide  (solid  line)  and  Model  Equation  Fit  (dashed  line).  Differences 
between  the  two  curves  are  highlighted  by  stippling  as  in  Figure  4 


where  u  is  the  windspeed.  Several  points  should  be  made  about  the  methodology 
described  in  Figure  6.  First,  the  Ri^  value  can  change  as  u  is  adjusted  at  each  grid 
point,  even  though  B  remains  constant.  Second,  the  use  of  Monin-Obukhov  scaling 
parameters  in  determining  n  is  consistent  with  using  the  same  scaling  relationships 
in  computing  B.  Since  the  assumption  of  a  constant  B  over  the  domain  is  only  an 
approximation,  we  allow  local  stability  effects  to  affect  the  model  wind  solution  using 
the  procedure  described  in  the  bottom  flowchart  of  Figure  6.  Thus  the  calculation 
of  local  grid  values  of  Ri,  ,  Q  ,  ,  and  n  allows  the  momentum  field  to  respond 

to  local  stability  changes  throughout  the  model  integration,  and  allows  changes  in 


the  momentum  field  to  feed  back  into  the  local  stability.  It  is  remarkable  that  the 
model  solution  of  Eq.  (15)  can  reach  a  minimum  constraint  within  10  to  20  steps 
despite  he  local  stability  adjustments  made  through  the  computations  described 
by  Figure  6. 
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Figure  6.  Flowchart  Describing  Calculation  of  Wind  Profile  Exponent  n  Based  on 
Values  of  Buoyancy  and  Windspeed.  Top  flowchart  shows  calculation  at  individual 
grid  point.  Bottom  flowchart  shows  how  exponent  calculation  leads  (o  changes  in 
momentum  flux,  and  changes  in  flux  feedback  to  alter  subsequent  computations 
of  n 


Since  the  wind  profile  exponents  computed  in  the  model  apply  to  conditions 

fairly  close  to  the  surface,  use  of  these  exponents  to  describe  wind  conditions  at 

higher  levels  in  the  surface  layer  leads  to  overestimates  of  the  wind  speed.  This 

occurs  because  the  value  of  n  in  Eq.  (21)  is  sensitive  to  choice  of  Z.  Field  observa- 
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tions  of  wind  profiles  documented  by  Gringorten  and  Grantham  and  many  others 
show  that  the  magnitude  of  the  wind  speed  increases  rapidly  over  the  lowest  10  m 
(see  Gringorten  and  Grantham's  Tables  6.  1,  6.  3,  and  6.  5).  The  use  of  model 
wind  profile  exponents  to  determine  windspeeds  at  higher  levels  is  further  dis¬ 
cussed  in  Section  5  and  Appendix  A,  and  it  will  be  shown  that  to  calculate  accurate 
surface  layer  wind  profiles,  the  calculation  height  Z  in  Eq.  (21)  must  be  changed 
to  represent  the  entire  surface  layer  more  accurately. 


2.3  Appropriateness  to  Mesoscole  Meteorological  Analysis 

One  question  that  surely  arises  from  the  mathematical  derivations  presented 
in  the  previous  two  sections  is:  Do  variational  principles  actually  describe  what 
the  lower  atmosphere  is  doing  over  complex  terrain?  The  answer  to  this  question 
lies  in  the  approach  taken  in  the  formulation  of  the  basic  model  equation  [Eq.  (15)]. 
Recall  that  the  equation  of  motion  can  be  derived  through  variational  techniques  by 
employing  several  simplifying  assumptions  such  as  incompressible  and  inviscid 
flow,  or  isentropic  motion  and  inviscid  flow.  These  approaches  all  retain  the 
essential  dynamics  as  embodied  in  the  equation  of  motion.  The  windflow  analysis 
is  based  on  kinematic  constraints,  and  represents  a  kinematic  analysis  of  steady - 
state  conditions  under  which  the  momentum  advection  and  buoyancy  force  approxi¬ 
mately  balance  under  the  constraints  of  topography  (through  the  terrain -following 
coordinates)  and  mass  conservation  (through  the  incompressibility  assumption). 
Therefore,  Eq.  (15)  essentially  describes  a  solution  of  atmospheric  flow  and  does 
not  attempt  to  describe  the  dynamic  processes  that  lead  to  the  solution.  In  many 
ways,  the  rmdel  solution  expressed  by  Eq.  (15)  is  similar  to  that  obtained  by  run¬ 
ning  a  PE  model  through  some  period  of  time  until  the  model  reaches  a  steady-state 
solution.  Several  modeling  studies  of  slope  flows  have  attempted  this,  and  reveal 
that  it  may  be  easier  to  attain  steady-state  conditions  under  stable  nocturnal  condi¬ 
tions  than  under  unstable  conditions.  During  the  night,  a  balance  can  occur  along 

slopes  between  local  cooling  due  to  radiational  and  turbulent  effects,  and  advective 
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warming  due  to  developing  stratification  (McNider  and  Pielke  ).  Even  so,  surges 


40.  Gringorten,  I.  I.  ,  and  Grantham,  D.  D.  (1983)  Winds  as  a  Function  of  Height, 
in  Winds;  Chapter  6,  1983  Revision,  Handbook  of  Geophysics  and  Space 
Environments,  AFGL-l'ft -83  -0080,  AD  A132018. 

4  1.  McNider,  R.T.,  and  Pielke,  R.A.  ( 1984)  Numerical  simulation  of  slope  and 
mountain  flows,  J.  Climate  and  Appl.  Meterol.  23;  144  1  - 1453. 
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in  downslope  flow  intensity  can  occur,  producing  local  accelerations  in  the  flow 

(see  the  analytical  study  of  McNider  ).  During  the  daytime  hours,  other  physical 

processes  such  as  horizontal  diffusion  of  momentum  by  turbulent  eddies  can 

account  for  accelerations  in  the  momentum  field  that  persist  for  several  hours 
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before  reaching  a  steady  state  (Banta  ). 

Several  modeling  studies  have  attempted  to  simulate  flows  over  complex 

terrain  using  variational  techniques.  Most  of  these  models  fall  into  the  so-called 

mass  conservation  type,  since  they  relay  on  the  continuity  equation  in  ore  form 

or  another  for  their  analysis.  An  early  attempt  at  this  type  of  model  was  made  bv 
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Anderson,  who  generated  terrain-induced  wind  fields  using  both  the  continuity 

equation  (assuming  incompressible  flow),  and  a  Poisson  equation  to  represent 

thermal  effects  such  as  the  urban  heat  island  and  topographic  maxima  or  minima. 

This  analysis  was  used  for  generating  regional  windflow  climatologies  for  air 
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pollution  applications.  Dickerson  developed  a  two-dimensional  mass -consistent 

2  g  <y  y 

analysis  based  on  a  variational  least  squares  technique  first  proposed  by  Sasaki, 

Dickerson's  analysis  sought  to  minimize  the  "error"  between  the  observations  of 

velocity  flux  components  and  height  of  the  inversion  base,  and  the  adjusted  values 

of  these  variables,  weighted  by  appropriate  factors  based  on  the  error  variance  of 

the  observed  field.  A  three-dimensional  version  of  this  model  was  also  developed 
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and  documented  by  Sherman.  Other  variational  analysis  models  were  proposed 
47  48 

by  Fox  et  al  and  Goodin  et  al,  for  windflow  modeling  problems,  and  also 
followed  Sasaki's  technique. 

All  of  the  above-mentioned  modeling  studies  attempted  to  produce  kinematically - 
constrained  windflow  analyses  by  assuming  some  type  of  steady-state  condition 

42.  McNider,  R.T.  (1982)  A  note  on  velocity  fluctuations  in  drainage  flows, 

J.  Atmos.  Sci.  39: 1658-  1660. 

43.  Banta,  R.  M.  (1986)  Daytime  boundary-layer  evolution  over  mountainous 

terrain.  Part  II:  Numerical  studies  of  upslope  flow  duration,  Mon.  Wea.  Rev. 
114:1112-1130.  - 

44.  Anderson,  G.  E.  (197  1)  Mesoscale  influences  on  windfield,  J.  Appl.  Meteorol. 

10:377  -386. 

45.  Dickerson,  M.  H.  (1978)  MASCON  -  A  mass-consistent  atmospheric  flux  model 

for  regions  with  complex  terrain,  J.  Appl.  Meteorol.  17 :24  1  -253. 

46.  Sherman,  C.  A.  (  1978)  A  mass -consistent  model  for  wind  fields  over  complex 

terrain,  J.  Appl.  Meteorol.  17  ;3  12 -3  19. 

47.  Fox,  D.G.,  Fosberg,  M.A.,  Marlatt,  W.  E.  ,  and  Reeser,  W.  (  1976)  Analysis 

of  mountain  air  quality,  Proc.  of  Third  Symposium  on  Atmospheric  Turbu¬ 
lence,  Diffusion,  and  Air  Quality,  Am.  Meteorol.  Soc.  ,  Boston,  Mass., 
pp.  470-475. 

48.  Goodin,  W.  R.  ,  McRae,  G..J.,  and  Seinfeld,  .1.11.  ( 1980)  An  objective  analysis 

technique  for  constructing  three-dimensional  urban-scale  windfields, 

J.  Appl.  Meteorol.  19:98-108. 


coupled  with  an  appropriate  condition  such  as  incompressibility  or  mass -consistent 
flow.  The  advantages  of  these  models  are  outlined  below: 

(1)  Simplicity  of  physics, 

(2)  Computational  efficiency  (they  can  be  run  on  smaller  machines 
than  PE  models), 

(3)  High  horizontal  and  vertical  resolution  over  small  regions, 

(4)  Inclusion  of  some  physical  processes  (though  kinematic  in  nature), 
thus  superior  to  objective  analyses,  and 

(5)  Applicability  to  data-sparse  areas  (related  to  No.  4). 

There  are  obvious  disadvantages  to  these  types  of  models.  Some  of  them  are  listed 
here,  and  are  further  described  throughout  this  report: 

(1)  Limited  applicability  to  only  certain  meteorological  situations, 
such  as  topographically -forced  flows, 

(2)  Non -predictive  nature, 

(3)  Emphasis  on  analysis  of  only  one  tvpe  of  variable  (usually  winds), 

(4)  Difficulties  in  determining  appropriate  means  of  initialization 
and  data  input,  and 

(5)  Limited  physical  paramete rizations. 

3.  HARDWARE  REQUIREMENTS  AND  MODELING 
SYSTEM  SOFTWARE 

This  section  describes  the  different  computers  on  which  we  run  the  windflow 
modeling  svstem.  We  also  outline  the  computer  architecture  for  the  windflow 
model  and  the  two  plotting  programs  used  to  display  the  output.  Complete  program 
and  subroutine  descriptions  are  included,  along  with  a  flowchart  of  the  model  pro¬ 
gram  . 

3.1  Hardware 

The  windflow  model  runs  on  the  Zenith  100  (Z-100)  and  Zenith  248  (Z-248) 
microcomputers,  and  on  the  Control  Data  Corporation  Cyber  860  main  frame 
com  puter. 

The  Z-100  on  which  most  of  the  prototype  development  was  done  is  a  16-bit 
desktop  machine  with  448  K  random  access  memory  (RAM),  an  added  co-processor 
("math  (  hip")  that  speeds  up  arithmetic  operations  by  about  10  times,  dual 
diskette  drives,  and  a  high -resoluf  ion  printer.  Screen  graphics  produced  by  the 
model  peripheral  routines  included  windflow  over  the  terrain  and  wind  profiles. 
These  were  output  to  the  12-in.  diagonal  monochrome  monitor  capable  of  handling 
640  '  22a  addressable  points  and  2.6  printable  lines  of  80  columns  each. 


The  Z-100  available  at  Ft.  Devens  did  not  have  a  math  chip.  Run  time  for  a 
15  X  15  array  (200-m  horizontal  grid  spacing)  on  this  Z-100  varied  from  12  to  20 
min,  depending  upon  the  atmospheric  stability  conditions. 

The  AFGL  Cyber  860  mainframe  computer  was  used  during  early  model  devel¬ 
opment  because  of  its  superior  size  and  compilation  speed.  The  Cyber  is  a  60-bit 
large-scale  machine  that  operates  on  a  24-hr  schedule.  As  a  result,  many 
model  runs  were  submitted  simultaneously  for  overnight  execution  during  the 
development  and  debug  stages  of  the  modeling  effort.  The  Cyber  computer  was 
also  used  for  developing  the  terrain  data  base  files  from  the  terrain  data  tapes. 

These  files  were  then  transferred  to  floppv  disks  for  use  on  the  Z-100  and  Z-248. 

During  model  development,  AFGL  acquired  a  Z-248  m icrocomputer.  The 
Z-248  is  IBM -compatible,  about  ten  times  faster  than  the  Z-100  (with  the  math  chip), 
has  expansion  memory  to  3200K  RAM  with  20  megabyte  hard  disk  storage,  a  360K 
diskette  drive,  and  superior  graphics  capabilities. 

3.2  Software 

The  software  for  the  windflow  model  was  coded  in  Zenith  Microsoft  Fortran  77. 
The  graphic  s  utilities  were  written  in  Microsoft  Z-Basic. 

3.2.1  WINDFLOW  MODEL 

The  main  program  MODEL  calls  three  major  modules  INPUT!),  WFMD  and 
FLOl'TD.  INPUT!)  calls  six  subroutines  in  the  process  of  setting  up  a  particular 
model  run.  These  input  such  variables  as  terrain  heights  for  the  selected  domain, 
date  and  time,  wind  speed  and  direction,  temperature,  and  cloud  cover.  Buoyancy 
and  stability  class  are  then  calculated.  WFMD  is  the  windflow  computational 
module,  calling  five  subroutines  to  produce  a  high  resolution,  two-dimensional 
analysis  of  surface  winds.  FLOUT!)  creates  three  output  files  that  contain  the 
wind  component  arrays  (u  and  v)  for  the  domain,  the  power-law  wind  profile 
exponents,  and  the  model  input  parameters. 

The  model  provides  spue  for  a  terrain  array  size  up  to  7  5  •  75  elements 
(44  ■  fi'1  on  the  Z-lOO).  The  model  requires  a  minimum  core  of  about  448K  to  be 
run.  The  size  is  located  in  a  PARAMETER  statement  in  an  "INCLUDE"  file  and 
can  be  altere  I  easily.  All  subroutines  access  common  variables  marie  available 
through  live  "INCLUDE"  files  pro\  ided  at  compilation  time.  A  list  of  those  files 
h 1 1  low  s : 

(1>  W  IN  It/./..  INC —contains  the  allocation  for  array  size  tor  the  model, 
used  lor  the  terrain  heights,  surface  roughness,  wind  components, 
wind  exponent,  and  buoyancy  arrays.  In  addition,  space  is  allocated 
for  the  computed  terrain  slopes,  and  common  variables  for  the  windflow 


(2)  WIN DBZ.  INC— provides  space  for  the  wind  speed,  direction  and 
surface  temperature. 

(3)  DIFUSAA.  INC —provides  space  for  constants  used  in  calculating 
buoyancy. 

(4)  WINDAZ.  INC— provides  common  storage  for  the  terrain  input  file  name. 

(5)  DATEZ.  INC —provides  space  for  date/ time  information. 


The  routines  called  by  the  modules  INPUTD,  WFMD,  and  FLOUTD  are  des¬ 
cribed  below.  A  complete  flowchart  of  the  windflow  model  is  displayed  in  Figure  7. 
INPUTD  executes  all  input  subroutines.  These  are: 


(1)  GETTER —inputs  domain  terrain  data  array  limits  (horizontal  and 
vertical),  rid  size,  longitude  and  latitude,  roughness  indicator. 

—  inputs  terrain  data 

—  creates  constant  roughness  array  (50  cm  over 
Ft.  Devens  in  prototype) 

—  sets  calculation  height  for  windflow  at  32.8  ft  (10  m)  AGL 

—  calculates  terrain  slopes 

—  inputs  date  and  time 

(2)  SOLAR— calculates  solar  radiation  based  on  date,  time,  latitude 
and  longitude. 

(3)  WIN  DIN  — inputs  wind  direction  in  degrees  and  speed  in  knots  and  uses 
these  to  initialize  wind  component  arrays  (u  and  v,  in  m/s)  at 

each  gridpoint.  Rotates  components  45’  to  orient  axes  along 
direction  of  calculated  terrain  slopes. 

—  inputs  cloud  cover  in  eights  (use  look-up  table  to  adjust 
solar  radiation  computation  for  type  of  cloud  cover) 

(4)  STAB3 —calculates  the  stability  parameter  (between  0.5-6.  0). 

(5)  TSRFIN  — inputs  surface  temperature  in  ’F  and  converts  it  to  Kelvin. 

(6)  BUOYIN —calculates  the  buoyancy  at  each  grid  point  using 
temperature  and  stability  parameter. 

—  displays  buoyancy  and  stability  parameter  and  returns 
to  main  program. 


WFMD  sets  the  maximum  number  of  relaxation  steps  at  60.  It  displays  ihe 
array  size  for  the  windflow  calculation  and  allows  choice  of  smaller  array  window. 

It  also  displays  the  current  relaxation  step  and  tests  to  see  if  constraint  has  reached 


minimum.  This  module  calls  five  subroutines. 
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Figure  7.  Flow  Chart  of  Windflow  Model  Program 
and  Peripheral  Software.  Also  shown  are  flowcharts 
for  the  main  modules  INPUTD  and  WFMD 


(1)  RELAX  — calls  routines  to  calculate  local  contributions  to  the  acceleration 
residual. 

—  accumulates  residuals 

—  calculates  partial  derivatives  of  residuals  with  respect  to 
wind  components  at  each  gridpoint 

—  applies  incremental  adjustments  to  wind  field 
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Figure  7.  Flow  Chart  of  Windflow  Model  Program  and  Peripheral  Software.  Also 
shown  are  flowcharts  for  the  main  modules  INPUTD  and  WFMD  (Contd) 


1  (2)  WINEXP  —  computes  the  wind  profile  exponent  at  each  gridpoint. 

(3)  SET  BOX— selects  the  parameters  for  the  local  residual. 

(4)  RESL.OC  — calculates  the  integral  of  the  residual  squared  over  a 
local  volume  element  (flux  box). 

|  (5)  OUTPUT —rerotates  calculated  wind  components  45°  to  return 

,  to  original  Cartesian  coordinates. 
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FLOUTD  creates  three  formatted  output  files  to  be  used  by  the  utility  plotting 
routines  to  display  output  on  the  screen. 

(1)  WIND.  OUT— contains  u,  v.  component  arrays  for  a  size  m  x  (west -east) 
by  m  y  (south -north). 

(2)  WEX.  OUT  —contains  wind  profile  exponent  array  for  m  x  by  m  y  size. 

(3)  MODIN. OUT—  contains  model  input  and  output  parameters.  These 
are  model  computation  height,  day,  month,  year,  time,  wind  direction 
and  speed,  surface  temperature,  cloud  cover,  cloud  type  if  cloud 
cover  is  7/8  or  8/8,  diagnosed  stability  parameter,  and  mean  wind 
profile  exponent  over  window. 

3.  2.  2  HORIZONTAL  WIND  PLOT  FOR  SCREEN 

This  routine,  called  WINDPLT,  displays  the  windfield  over  a  selected  window 
(the  user  can  choose  the  array  size  for  the  display).  The  winds  can  be  displayed 
either  as  arrows  with  length  proportional  to  speed  (in  this  option  the  program  cal¬ 
culates  the  maximum  windspeed  encountered  over  the  chosen  display  window  and 
asks  the  user  if  he  would  like  to  use  this  value  for  the  scale),  or  in  conventional 
plot  notation  with  barbs.  The  program  displays  the  winds  in  kt,  and  also  displays 
the  date/time,  mean  direction  and  speed  throughout  the  window,  mean  profile 
exponent,  stability  parameter,  and  model  calculation  height  in  feet.  Depressing 
the  space  bar  on  the  keyboard  after  the  display  is  completed  causes  the  program 
to  ask  if  the  user  wants  another  plot  or  to  end  the  program.  There  are  two  ver¬ 
sions  of  this  program;  one  for  the  Z-100  and  another  for  the  Z-248. 

3.  2.  3  VERTICAL  WINDSPEED  PROFILE  SCREEN  PLOT 

This  routine,  called  PROFILE,  displays  a  windspeed  profile  from  32.8  ft  to 
300  ft  if  the  user-  selects  the  speed  in  kt  option,  or  from  10  m  to  100  m  if  the  ms  ' 
option  is  selected.  The  user  inputs  the  windspeed  and  exponent  using  either  the 
model -generated  values  obtained  from  the  horizontal  wind  plot,  or  user -selected 
values.  As  in  the  horizontal  plot,  hitting  the  keyboard  space  bar  results  in  a 
choice  of  another  plot  or  program  termination.  This  program  also  exists  in  Z-100 
and  Z-248  versions. 

4.  A  USER'S  GUIDE  TO  THE  WINDFLOW  MODELING  SYSTEM 

The  purpose  of  this  chapter  is  to  familiarize  the  user  with  the  operation  of  the 
windfiow  model.  Since  a  large  group  of  potential  model  users  mav  be  operational 
forecasters  with  little  training  in  boundary -laver  meteorology,  we  have  included 
Section  4.  1  to  give  the  forecaster  an  introduction  to  the  dynamics  of  the 


lower  atmosphere.  This  introduction  will  serve  as  a  basis  for  subsequent  dis¬ 
cussions  of  the  model's  formulation  (Section  4.2)  and  rules  for  interpreting  the 
output  (Section  4.3)  under  varying  atmospheric  conditions.  Becoming  familiar 
with  this  section  will  make  the  interpretation  of  the  results  of  the  Ft.  Devens  field 
tests  in  Section  5  easier. 

4.1  Introduction  to  the  Dynamics  of  the  Lower  Atmosphere 

This  section  is  not  intended  to  provide  a  comprehensive  review  of  the  dynamic 
meteorology  of  the  lower  atmosphere.  Instead,  it  provides  the  necessary  back¬ 
ground  for  understanding  how  the  windflow  model  operates,  as  well  as  the  limita- 
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tions  it  has.  Material  in  this  section  is  drawn  primarily  from  Anthes  et  al  and 
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Wyngaard. 

The  windflow  model  produces  a  high-resolution  analysis  of  surface  winds  as 
influenced  by  the  stability  of  the  lower  atmosphere  and  the  topography.  To  better 
understand  the  relationship  described  by  the  model,  we  begin  with  some  definitions. 
The  lower  atmosphere  is  often  referred  to  as  the  planetary  boundary  layer,  or 
PBL.  The  PBL  is  defined  as  that  layer  of  the  atmosphere  where  the  effects  of  the 
earth's  surface  are  dominant.  During  the  day,  the  PBL  (typically  1-2  km  deep) 
is  influenced  by  turbulent  mixing  of  momentum,  heat,  and  moisture.  This  turbu¬ 
lent  mixing  occurs  because  the  atmosphere  is  largely  transparent  to  solar  radia¬ 
tion  and  most  of  the  sunlight  heats  the  ground  instead  of  the  air.  As  the  ground 
heats  up  during  the  day,  it  conducts  heat  both  downward  and  upward.  The  upward 
conduction  causes  the  air  very  close  to  the  ground  to  heat  up  and  expand,  thus 
rising  above  cooler,  denser  air  above  it.  Bubbles  of  heated  air  begin  to  rise  while 
cooler  air  sinks  to  the  ground  and  is  heated  by  the  ground.  As  these  heated  bubbles 
of  air  rise,  they  mix  with  air  aloft  and  heat  it.  The  rising  and  sinking  bubbles  are 
called  convection  currents,  and  they  are  the  primary  means  by  which  turbulent 
mixing  occurs  in  the  daytime  PBL.  At  night,  the  ground  radiates  heat  away  to 
space,  and  the  air  above  the  ground  slowly  cools  by  radiation  too.  The  PBL  be¬ 
comes  much  shallower  because  turbulent  mixing  due  to  surface  heating  has  ceased 
and  a  surface -based  temperature  inversion  forms.  If  there  is  mixing  at  night,  it  is 
usuallv  caused  by  strong  winds  (mechanical  mixing).  The  depth  of  the  PBL  in¬ 
creases  with  windspeed  and  surface  roughness. 

49.  Anthes,  R.  A.  ,  Panofskv,  H.  A.  ,  Cahir,  J.  .1.  ,  and  Rango,  A.  (  1978)  The 

Atmosphere  (2nd  Ed.),  Merrill  Publishing  Co.  ,  Columbus,  Ohio,  pp.  128-131. 

50.  Wvngaard,  .J.C.  (1985)  Structure  of  the  planetary  boundary  layer  and  implica¬ 

tions  for  its  modeling,  .J.  Climate  and  Appl.  Meteorol.  24:1  13  1  -1  142. 


Typical  profiles  of  temperature  for  daytime  and  nighttime  conditions  are 
shown  in  Figure  8.  Notice  that  the  effect  of  strong  winds  is  to  alter  the  thermal 
structure  through  mechanical  mixing  so  that  the  daytime  and  nighttime  profiles 
become  more  similar. 

The  windflow  model  considers  that  portion  of  the  PBL  known  as  the  surface 
laver.  The  surface  layer  is  usually  considered  to  be  the  lowest  10  percent  of  the 
PBL.  It  is  a  layer  through  which  vertical  transport  (flux)  of  momentum,  heat,  and 
moisture  can  be  considered  constant.  During  the  daytime  the  surface  layer,  like 
the  PBL,  is  well-defined.  The  problems  associated  with  defining  the  nighttime  PBL 
also  applv  to  the  surface  layer  (this  problem  will  be  highlighted  further  in  Section 
-4.  3).  The  temperature  profiles  over  the  lowest  100  m  in  Figure  8  are  typical  of 
the  surface  laver.  Notice  that  strong  vertical  temperature  gradients  can  exist  in 
the  surface  laver.  The  surface-layer  stability  is  a  measure  of  the  magnitude  of 
the  vertical  temperature  gradient.  Stability  used  in  this  context  is  not  the  same 
as  the  "static  stability"  used  to  describe  conditions  favorable  for  thunderstorm 
formation.  There  are  three  basic  states  of  surface -layer  stability  based  on  the 
vertical  temperature  gradient.  Unstable  conditions  (that  usually  appear  during 
the  dav)  occur  when  the  vertical  temperature  gradient  is  less  than  -1°C  100  m  *. 
Stable-  conditions  (usually  appearing  at  night)  occur  when  the  vertical  temperature 
gradient  is  greater  than  -PC  100  m  \  Neutral  conditions,  that  usually  occur 
under  cloudy  and/or  windy  conditions,  or  at  sunrise  and  sunset,  correspond  to  a 
vertical  temperature  gradient  equal  to  -1°C  100  m  The  -1°C  100  m  1  value  is 
also  known  as  the  drv  adiabatic  lapse  rate  (the  temperature  lapse  rate  obtained 
bv  following  the  drv  adiabats  on  a  conventional  Skew  T -Log  P  diagram).  It  is 
important  to  understand  these  definitions  concerning  surface -layer  stability  since 
thev  are  not  normally  resolved  on  conventional  temperature-dewpoint  soundings 
(for  instance,  an  unstable  layer  on  a  sounding  appears  as  a  "superadiabatic"  layer, 
which  does  not  occur  very  often  since  the  vertical  resolution  is  not  fine  enough). 

For  this  reason,  the  model  diagnoses  stability  based  on  the  surface  input  informa- 


Figure  8.  Profiles  of  PBL  Temperature  for  Daytime  and  Nighttime  Conditions,  Under  Light  Wind 
Clear  Sky  Conditions  (solid  lines)  and  Cloudy  Skies/Strong  Winds  (dashed  lines) 


Typical  surface -layer  wind  profiles  for  daytime  and  nighttime  conditions  are 
shown  in  Figure  9.  During  the  daytime,  windspeed  increases  with  height  to  about 
100  m  and  is  approximately  constant  throughout  the  rest  of  the  PBL.  At  night, 
the  winds  are  typically  light  near  the  surface,  but  increase  rapidly  with  height 
until  the  top  of  the  nocturnal  temperature  inversion  is  reached.  At  this  level, 
a  wind  maximum  called  a  "nocturnal  jet"  is  often  encountered.  Quite  often,  the 
nocturnal  jet  is  referred  to  as  a  "low-level  jet.  "  These  two  phenomena  are  closely 
related,  but  the  low-level  jet  is  a  broader  category  that  includes  those  jets 


whose  formation  can  be  explained  by  their  interactions  with  upper-level  jets 

5  1 

(Uecellini  and  Johnson  )  as  well  as  those  strong  low-level  jets  associated  with 
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nocturnal  thunderstorm  activity  over  the  Great  Plains  (Means;  Lettau; 
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Blackadar;  Pitchford  and  London;  Bonner  and  many  others).  The  nocturnal 


jet,  bv  contrast,  may  not  be  that  strong,  and  is  usually  associated  with  the  forma¬ 
tion  of  the  nocturnal  temperature  inversion.  It  should  be  noted  that  under  neutral 
conditions  both  the  wind  and  temperature  profiles  tend  to  be  more  uniform  with 
height  due  to  the  effects  of  mechanical  mixing.  Now  that  we  have  discussed  the  PBL 
and  surface  layer,  we  can  summarize  the  wind  and  temperature  characteristics  for 
each  of  the  three  stability  classes.  This  is  done  in  Table  2. 

Wo  now  discuss  the  effects  of  topography  on  the  surface  windflow.  In  general, 
the  profile  relationships  outlined  in  Table  2  are  still  valid  over  sloping  terrain, 
but  with  some  important  additional  effects. 
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ndspeeds  for  Daytime  and  Nighttime  Conditions 


Table  2.  Wind  and  Temperature  Characteristics  Under  Different  Stability 
Conditions 


Surface -Layer 
Stability 

Most  Common 
Meteorological 
Conditions 

Wind  Profile 
Characteristics 

Temperature 

Profile 

Characteristics 

Unstable 

Daytime*,  Light  Winds, 
Clear  to  Scattered 

Skies 

Increase  w/height 
in  lowest  50-100  m, 
approximately 
uniform  to  top 
of  PBL 

Strong  decrease 
w/height  in 
lowest  50-100  m, 
then  dry  adiabatic 

(-1°C  100  m"1) 
above  to  inversion 
at  top  of  PBL 

Stable 

Nighttime**,  Light 
Winds,  Clear  to 
Scattered  Skies 

Near  calm  at 
surface,  strong 
increase  w/height 
above  surface 
to  top  of 
temperature 
inversion,  then 
decrease  above  to 
daytime  top  of 

PBL 

Increase  or 
isothermal 
w/height  to  a 
level  somewhere 
below  daytime 
top  of  PBL, 
then  dry  adiabatic 
above  to 
inversion  at 
daytime  top  of 

PBL 

Neutral 

Sunrise,  Sunset, 

Strong  Winds  and/or 
Cloudy  Conditions 

Increase  w/height 
in  lowest 

50-100  m,  then 
uniform  above  to 
top  of  PBL 
(daytime)  or  to 
stable  layer 
(nighttime) 

Approximate 
dry  adiabatic 
w/height 
< - 1  °C  100  m'1 ) 
to  top  of  PBL 
(daytime)  or 
to  stable  layer 
(nighttime) 

An  unstable  layer  can  also  be  found  near  the  surface  in  urban  areas  at  night 
(known  as  urban  heat  island  effect). 


A  stable  layer  near  the  surface  can  sometimes  be  found  during  the  day  in 
winter  at  high-latitude  locations,  or  at  mid-latitude  locations  with  either 
a  snow  cover  or  arctic  high  conditions. 


During  the  day,  the  ground  is  heated,  but  the  heating  is  stronger  over  sloping 
terrain  facing  the  sun  than  over  level  ground  due  to  a  more  direct  solar  angle.  The 
slope  surface  heating  induces  a  horizontal  temperature  gradient  between  the  air 
close  to  the  slope  and  the  free  atmosphere  at  the  same  elevation  (see  Figure  10). 
These  temperature  differences  induce  differences  in  pressure,  so  a  pressure 
gradient  force  from  the  free  air  to  the  slope  is  generated.  Meanwhile,  an  unstable 
vertical  temperature  gradient  exists  above  the  sloping  ground  and  produces  an 


upward  buoyancy  force.  These  two  forces  cause  the  air  to  accelerate  horizontally 
toward  the  slope  and  vertically  above  the  slope,  and  the  presence  of  the  terrain 
serves  to  direct  this  accelerated  flow  "upslope."  At  night,  the  pressure  gradient 
force  is  directed  away  from  the  slope  (since  the  air  temperature  above  the  sloping 
ground  is  cooler  than  the  free  air  at  the  same  level),  and  the  stable  lapse  rate 
above  the  sloping  ground  causes  a  downward  buoyancy  force.  The  combination  of 
these  two  forces  in  the  presence  of  the  terrain  causes  a  "downslope"  wind  to  appear. 


DAYTIME 


RESULTING  UPSLOPE  FLOW' 


SLOPE  ANGLE - 


T  HEIGHT 

Ta  >Tb,  SO  Pa<Pb 

^  T  =  TEMPERATURE 
P  =  PRESSURE 

PGF  =  PRESSURE  GRADIENT 
FORCE 

BF  =  BUOYANCY  FORCE 


NIGHTIME 


RESULTING  DOWNSLOPE  "T  ^ 
FLOW  f  BF 


HEIGHT 


ta <  Tb,  so  pa>pb 


SLOPE  ANGLE - 


Figure  10.  Diagrams  Showing  Slope  Heating  (top)  and  Cooling  (bottom)  Effects 
on  Formation  of  Slope  Winds  (after  Atkinson^O) 


The  slope  flows  described  above  are  also  influenced  by  the  strength  of  the  pre¬ 
vailing  winds  as  well  as  the  presence  of  clouds.  Under  strong  winds  and/or  cloudy 
conditions,  the  slope  heating  effects  are  weakened  and  the  flows  along  the  slopes 
are  more  aerodynamic  (that  is,  influenced  by  the  forces  resulting  from  the  wind- 
flow  interaction  with  the  terrain  obstruction). 


4.2  Model  Formulation  and  Physical  Assumptions 

The  windl'low  model  runs  on  a  grid  system  that  produces  calculations  of  wind 
speed  and  direction,  taking  the  slope  of  the  surrounding  terrain  into  account.  The 
model  begins  with  a  uniform  windfield,  and  makes  small  adjustments  at  each  grid 
point  over  the  model  domain.  This  adjustment  is  called  a  relaxation  step.  After  a 
number  of  relaxation  steps,  the  model  has  integrated  the  effects  of  ambient  flow, 
stability,  and  topography  to  reach  a  solution  and  the  model  run  ends.  The  model 
may  take  a  different  number  of  steps  to  reach  a  solution  depending  upon  the  meteor¬ 
ological  situations.  At  the  end  of  this  section,  we  describe  a  way  for  the  user  to 
tell  approximately  how  many  steps  a  model  run  will  take  for  different  cases. 

The  easiest  way  to  describe  the  model  formulation  is  to  outline  a  sample 
session  with  the  model  (see  Figure  11).  The  Z-100  and  Z-248  versions  are  essen¬ 
tially  identical,  so  the  session  described  here  is  applicable  to  either  machine. 

The  model  first  reads  in  digital  terrain  elevation  data  for  the  area  of  interest,  noting 
the  latitude  and  longitude  of  the  location  for  calculating  solar  elevation  angles  and 
solar  radiation.  The  model  also  reads  in  the  surface  roughness,  a  parameter  that 
describes  the  character  of  the  land  surface.  The  model  calculation  height  is  set 
equal  to  10  m  above  ground  level  (AGL)  because  the  flow  at  this  level  is  strongly 
influenced  by  the  land  surface  features.  Next,  the  date  and  time  (in  GMT)  are  input 
for  the  model  to  complete  its  solar  radiation  computation.  The  wind  direction  and 
speed  are  entered  next  and  the  model  sets  up  the  uniform  initial  windfield.  If  light 
and  variable  or  calm  winds  are  reported  or  forecast,  use  either  the  predominant 
observed  direction  or  make  multiple  runs  using  different  directions  (the  default 
windspeed  in  the  model  is  1  kt).  The  next  piece  of  information  is  the  cloud  cover 
in  eighths.  This  input  presents  a  problem  when  cloud  decks  exist  at  different  levels. 
Table  3  describes  the  cloud  input  using  observations  and  military  terminal  aero¬ 
drome  forecast  (TAFs).  For  7/8  and  8/8  cloud  cover  (overcast  conditions),  the 
model  asks  you  ;  •>  choose  from  eight  different  cloud  types,  ranging  from  cirrus  to 
fog.  Table  4  shows  the  guidelines  for  input  of  cloud-type  based  on  either  observa¬ 
tions  or  TAF's.  Notice  that  for  broken  cloud  decks  with  an  overcast  layer  above, 
the  dominant  cloud  type  corresponds  to  the  deck  at  ceiling  level,  not  overcast  level. 
This  is  because  the  model  considers  the  effects  of  all  cloud  layers  on  the  amount  of 
solar  radiation  reaching  the  surface.  At  night,  the  downward  heat  flux  reaching  the 
surface  is  more  affected  by  the  presence  of  low  clouds  than  by  upper-level  cloud 
layers.  The  model  finally  asks  for  surface  temperature,  and  then  displays  the 
stability  parameter  corresponding  to  the  input  data  given  to  it  by  the  user.  In 
general,  a  stability  parameter  below  3.3  corresponds  to  unstable  conditions,  from 
3.  3  to  3.  7  is  neutral,  and  above  3.7  is  stable.  One  has  the  option  of  running  the 
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Figure  11.  A  Sample  Interactive  Session  Using  Surface  Observation  as  Input  to 
Wind  flow  Model.  Input  bv  user  is  indicated  by  an  "i"  in  left-hand  column  with 
a  "CR"  after  input  data  to  indicate  carriage  return  on  keyboard 


Table  3.  Cloud  Information  Input 


Using  Observations 

Observed  Sky  Condition 

Cover  Input  to  Model 

Clear 

0 

Scattered  (one  deck) 

1,2 

Scattered  (two  decks) 

3 

Scattered  (more  than  two  decks, 
no  ceiling 

4 

Broken  (ceiling  above  10,  000  ft, 
no  overcast  deck) 

5 

Broken  (ceiling  below  10.  000  ft, 
no  overcast  deck) 

6 

Broken  with  overcast  above,  or 

Fog  with  visibility  2  1  mi, 

<  3  mi* 

7 

Overcast,  or  Fog  with  visibility 
below  1  mi* 

8 

Using  TAP's: 

Input  cloud  cover  in  eights  from  sky  condition  on  TAF  using 
guidelines  above  for  multiple  cloud  decks. 

*For  situations  with  both  clouds  and  fog,  run  the  model  twice, 
using  the  cloud  condition  and  then  using  the  fog  condition. 


model  for  the  entire  domain  or  over  a  window  of  that  domain.  The  example  in  Fig¬ 
ure  11  shows  that  a  43  X  68  array  has  been  specified  for  Ft.  Devens,  Massachusetts, 
that  translates  into  2924  grid  points  at  which  model  calculations  are  performed. 

A  window  of  the  domain  has  been  chosen  that  corresponds  to  the  drop  zone  at 
Ft.  Devens.  This  window  is  a  15  X  15  array,  or  225  grid  points. 

The  number  of  steps  the  model  will  take  to  reach  a  final  solution  is  primarily 
dependent  upon  two  factors.  The  first  factor  is  the  stability  parameter,  and  the 
second  one  is  the  ruggedness  of  the  terrain.  The  model  will  converge  '  o  a  solution 
faster  if  the  stability  parameter  is  close  to  neutral  (between  3.3  and  3.7).  As  the 
stability  parameter  increases  bevond  3.7,  or  decreases  below  3.3,  the  model  run 
time  will  increase  as  well.  For  example,  over  Ft.  Devens,  the  model  took  about 
18  to  20  steps  to  reach  a  solution  for  very  stable  cases  (stability  parameter  -  6.0, 


Table  4.  Cloud  Type*  Input  (7 / 8  or  8/8  Cloud  Cover) 


Observations: 

The  predominant  cloud  deck  should  be  the  one  at  which 
the  ceiling  is  located  (that  is,  35  BKN  80  OVC  should  be 
treated  as  7/8  cloud  cover  with  predominant  cloud 
type  Sc). 

TAFs: 

Same  as  observations  (that  is,  2  St  008  5  Sc  035  7  AS 
100  CIG  035  should  be  treated  as  7/8  cloud  cover  with 
predominant  cloud  type  Sc) 

$ 

See  Appendix  B  for  definitions  of  the  cloud  type 
abbreviations. 


the  model  usually  runs  in  a  minimum  of  two  steps.  Over  more  rugged  terrain  such 
as  Vandenberg  AFB,  California,  the  model  run  time  increases  for  stabilities  just 
below  3.3  and  above  3.7,  but  remains  about  the  same  for  the  extreme  cases.  This 
relationship  is  described  in  Figure  12. 

4.3  Rules-of-Thumb  for  Interpreting  Model  Output 

In  this  section,  we  document  those  situations  that  the  windflow  model  handles 
well,  as  well  as  those  cases  in  which  the  model  output  may  not  be  as  reliable.  The 
windflow  model  output  can  be  displayed  in  two  wavs.  The  first  is  bv  a  horizontal 
plot  of  the  winds  at  10-m  AGL.  It  is  advisable  to  overlay  a  terrain  map  with  geo¬ 
graphic  features  on  the  horizontal  wind  plot  in  order  to  see  the  terrain -induced  wind 
features  in  graphic  detail.  The  second  is  a  vertical  windspeed  profile  over  the 
lowest  100  m  (about  300  ft),  described  in  Section  3.3.  This  wind  profile  output 
warrants  separate  discussion  that  appears  in  Section  4.3.2. 

4.3.  1  INTERPRETATION  OF  HORIZONTAL 
WIND  PLOTS 

In  general,  the  windflow  model  produces  the  best  results  when  it  is  initialized 
with  input  information  that  is  representative  of  the  area  being  analyzed.  The  model 
domain  is  usually  on  the  order  of  about  10  *'  10  km;  so  a  single  observation  is  some¬ 
what  representative.  However,  we  should  give  vou  a  few  pointers  to  ensure  that 
this  is  true.  First,  a  good  knowledge  of  Ihe  microclimate  is  essential  before 
running  the  model.  Consider  such  factors  as  proximity  of  water  bodies  such  as 


MODEL  RELATION  STEP  AS  A  FUNCTION  OF  TERRAIN  RUGGEDNESS 

AND  STABILITY  PARAMETER 


UNSTABLE  NEUTRAL  STABLE 

STABILITY  PARAMETER  — 

Figure  12.  Diagram  Illustrating  Approximate  Number  of  Steps  the  Model  Takes 
tn  Reac  h  a  Solution  as  a  Function  of  Stability  Parameter  and  Terrain  Ruggedness. 

I  hi.-  diagram  should  be  used  as  a  guide  only;  actual  results  could  be  different 

large  lakes  or  oceans  that  could  produce  local  circulations  such  as  sea  and  land 
breezes.  These  loc  al  phenomena  may  produce  different  stabilities  in  different 
locations  (for  example,  neutral  conditions  at  the  shoreline  and  unstable  inland),  and 
using  the  observed  wind  from  one  locale  (sut  h  as  the  beach)  may  not  be  representa- 
1 1 v e  of  the  entire  area  vou  are  interested  in.  Multiple  model  runs  should  be  made 
with  different  tvpe  of  inputs  that  cover  the  different  possibilities  for  a  given  time. 
1'smg  the  coastal  example,  try  using  input  conditions  (such  as  neutral  for  the 
i  oastline,  then  unstable  for  the  interior),  that  will  yield  different  model  solutions 
over  the  domain,  and  compare  the  results.  If  vou  are  unsure  of  the  representa¬ 
tiveness  of  the  local  wind  observation,  compare  it  against  a  derived  value  from  the 
local  sin  face  analysis,  such  as  a  geostrophic  or  gradient  wind.  Make  model  runs 
using  these  different  wind  inputs  and  compare  the  results.  For  forecast  applica¬ 
tions,  the  output  Irom  a  model-output  statistics  package  can  give  you  a  representa¬ 
tive  ’observation"  to  use  in  the  model  run.  Compare  this  output  with  that  using  the 
TA  F  as  input. 
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II  i\  .  r.g  I  :!><•  different  wavs  in  which  to  arrive  at  a  representative 

.•  •  •:  *  I  a-  ::  ->  i<  1,  At-  >  I  i  mil  bsi  uss  how  to  interpret  the  output.  If  multiple 
it  .  ";n-  t.’f  ::..t  a-  rm  a  *»  i\  en  situation,  the  user  must  choose  which  output  is 
•  :  .  Ora  A  t.  lo  tins  is  to  have  alreadv  prepared  a  model  wind  climatology 

-.’..it)’:-  pa  ferns  at  different  times  of  the  dav,  times  of  the  year,  and 
A"-  t  run.  me  ‘he  m  >drl  f.»t  a  particular  forecast  time,  compare  the  output 
" f  'no-  •  h.  .  limand. >121,  al  simulation  for  that  tvpe  of  synoptic  flow  regime, 
or  :i  ii  .-o',  -h  >ul  i  gi\e  vt >u  a  first  guess  at  how  representative  the  model  run 
.  !>*  •  h*  nar’ii  ular  date-  and  time  of  interest. 

1  !.«•!•••  ar>-  t w  >  other  factors  necessarv  to  interpret  the  model  output.  One 
e  >  -  .'"in  Knowledge  ot  the  local  climate,  as  previously  discussed  with  regard 
he  :  e;>ri  -.-lit at  iveness  >t  the  input  observation.  The  second  factor  is  a  know- 
.'  o'  -he  >r,herent  -trengths  and  weaknesses  of  the  windflow  model.  Knowledge 
>  t*'l  '’lengths  and  weaknesses  is  important  regardless  of  the  location  and 
at  i>  >!:.  ani  potential  users  should  become  very  familiar  with  this  information 
>re  n.'itig  ’he  model.  Table  5  lists  the  model  strengths  and  weaknesses,  and 
now  iisi  uss  these  in  detail. 


T  able  5.  Model  Strengths  and  Weaknesses 


Strength 


—  Cold-air  drainage  flows 

—  Davtime  upslope  flows* 

—  Flow  under  homogeneous  wind  or  stability  conditions 

—  Effects  of  vegetation  on  the  flow  (light  winds,  stable  or 

unstable  conditions) 

Weaknesses 

—  Flows  under-  nonhomogeneous  wind  and  stability  conditions 

—  l-'lows  over  flat  terrain  areas 

—  Light  wind  cases  when  no  vegetation  data  are  available 

—  Model  cannot  be  applied  over  very  mountainous  regions  (see  text) 


There  is  a  tendency  for  the  model  to  exaggerate  upslope  flow 
conditions  over  gentle  topography  under  very  unstable  (for 
example,  stability  parameters  0.  5  -  1.0)  conditions 
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The  model  can  simulate  terrain -induced  flows  such  as  cold-air  drainage  and 
upslope  flows,  and  integrates  the  effects  of  stronger  winds  into  the  solution.  Strong 
flows  tend  to  show  little  terrain  influence  over  gentle  topography,  but  some  influence 
almost  alwavs  remains  over  rugged  terrain.  If  stability  conditions  are  fairly 
homogeneous  over  the  area,  the  model  results  can  be  relied  upon.  If  digital  vege¬ 
tation  height  data  are  available  for  a  location,  the  model  is  especially  capable  of 

integrating  the  effects  of  vegetative  cover  into  the  How  patterns  under  light  winds 
■>  3 

(see  I.anicci  ).  It  the  model  is  used  to  support  NBC  operations  (see  Section  1.  1), 
we  strongly  recommend  the  acquisition  of  vegetation  data  for  addition  to  the  terrain 
information  alreadv  used  bv  the  model. 

The  model  has  several  weaknesses  that  must  be  considered  when  interpreting 
the  uitpiit.  l'he  first  of  these  weaknesses  pertains  to  the  presence  of  non- 
homogeneous  stability  and  wind  conditions  over  an  area.  Several  examples  are 
g  iv  en  below : 

(1)  Coastal  regions  (especially  mountainous  areas  like 

\  andenberg)  —  part  icularlv  around  sunrise  and  sunset,  when  coast¬ 
line  stability  mav  be  neutral  or  slightly  stable,  while  the  interior 
remains  unstable;  sea  breezes  with  weak  gradient -level  winds,  or 
when  a  sea-breeze  "front"  produces  mesoscale  convergence  and 
localized  clouds  or  precipitation. 

(2)  \rrv  mountainous  regions  (areas  like  the  Continental  Divide  of 
North  America)  — In  these  locations,  the  high  terrain  is  under  a 
vastly  different  flow  regime  than  the  lowlands  are,  so  a  three- 
dimensional  model  would  be  more  appropriate. 

Idle  model  has  problems  over  flat  terrain  (although  most  regions  of  application  will 

have  few  flat  areas  in  them)  bec  ause  the  model  also  adjusts  the  winds  here,  leading 

to  directional  biases  in  the  output.  This  problem  is  not  a  serious  one  because  the 

model  run  usually  ends  before  the  bias  begins  to  appear  in  the  windfield  (Lanicci 
24  23 

and  Weber  ).  I.anicci  showed  that  the  model  produces  different  windfields  for 
vegetated  and  nonv  egeta  ted  areas  under  light  wind  conditions.  If  no  vegetation 
data  are  available  for  a  location,  the  model  wind  output  can  be  suspect  for  light 
wind  cases . 

4.  3.2  INTERPRETATION  OK  VERTICAL  WINDSPEED 
PROFILES 

Idle  vert  i,  al  w  mdspeed  profile  produced  by  the  model  was  briefly  discussed 
m  Sc.  turn  3.2.3.  Recall  the  surface -layer  wind  profiles  shown  in  Figure  9.  The 
model  profiles  are  produced  over  the  surface  layer,  which  is  assumed  to  be  100  m 
deep  luring  the  daytime  (on  the  average).  During  the  day,  the  windspeed  increases 


rapidlv  with  height  over  the  lowest  30  to  40  m,  then  begins  to  level  off  from 
40  to  100  m.  Since  the  com putat ion  height  is  equal  to  10  m,  the  model -derived 
profile  reflects  this  rapid  increase  with  height,  and  maintains  it  through  the  entire 
100 -m  laver  (see  example  from  model  output,  shown  in  Figure  13).  This  model¬ 
generated  profile  has  a  tendency  to  overestimate  the  winds  at  levels  above  10  m, 
as  will  be  documented  in  Section  5.  An  alternate  method  of  generating  these  pro¬ 
files  based  on  a  computation  height  of  31.  6  m  (the  geometric  mean  of  10  and  100  m, 
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recommended  bv  Panofsky  and  Dutton  ),  will  be  discussed  in  Section  5.  Those 
users  who  have  a  need  for  vertical  windspeed  profiles  should  read  the  next  section. 

At  night,  it  is  not  a  straightforward  task  to  produce  representative  wind 
profiles.  Recalling  Figure  0  again,  the  vertical  windspeed  profile  is  very  sharp, 
and  depends  upon  the  level  of  the  nocturnal  jet  (for  which  the  model  cannot  account). 
Therefore,  except  for  neutral  conditions  (for  which  the  winds  are  fairly  uniform 
with  height  even  at  night),  the  model  wind  profiles  are  very  suspect  at  night.  An 
example  of  such  a  nighttime  profile  generated  bv  the  model  is  shown  in  Figure  14. 
Notice  that  even  for  light  surface  winds  of  2-3  ms  \  the  model  profile  product's 
win  (speeds  in  excess  of  20  ms  '  only  100  m  above  the  surface.  We  present  the 
following  guidelines  for  use  of  the  model  profiles  under  nighttime  (stable)  conditions 

(It  Check  the  winds  and  temperatures  on  the  latest  available  sounding  to 

get  estimates  of  low-level  winds  and  location  of  the  temperature  inversion. 
Also  check  other  data  sources  such  as  pilot  balloons,  pilot  reports, 
numerical  model  products,  and  so  on,  for  low-level  winds. 

12)  In  votir  local  area,  watch  for  signs  of  the  presence  of  a  temperature 
inversion  (light  surface  winds,  clear  skies,  and  so  on).  During  the 
colder  months,  a  snow  cover  will  help  the  inversion  formation 
through  strong  radiational  cooling  after  sunset.  Watch  for  other  signs 
such  as  smoke  from  fireplaces  and  woodstoves  beginning  to  reduce 
visibility.  Also,  smoke  from  tall  stacks  w  ill  show  a  pattern  where  it 
rises  and  tans  out  when  inversions  arc  present.  In  the  carlv  morning, 
the  presence  >!  a  shallow  haze 'smoke  laver  is  another  good  visual 
mdii  ator  ol  die  inversion  laver. 

(3l  Alter  perlorming  steps  1  and  2.  run  the  mode!  and  plot  the  profile.  II 
die  inversion  height  is  near  100  m.  the  trend  ol  the  model  wind  profile 
tnav  be  reliable,  and  perhaps  just  the  spec  t  at  the  top  of  the  prolile 
needs  '■>  be  ie  lin  e  I.  II  the  top  of  the  inversion  laver  is  more  like 
200  m  or  greater,  the  prolile  is  probable  inaccurate,  especially  it 
the  vs  in  Is  pee  d  at  the  top  ol  the  prolile  exceeds  the  w  in  Ispeeds 
obt  a  met  in  step  1  . 
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5.  APPLICATION  OF  THE  MODELING  SYSTEM  TO  WEATHER 

SUPPORT  OPERATIONS  AT  FT.  DEVENS.  MASSACHUSETTS 

This  section  documents  the  field  test  we  conducted  at  the  Air  Weather  Service 
(AWS)  Detachment  at  Ft.  Devens,  Massachusetts  from  mid-September  to  early 
November  1986.  We  first  discuss  the  organization  and  operational  mission  of 
Detachment  12  (Del  12),  then  the  forecaster  training  that  was  conducted.  We 
next  describe  die  model  windflow  climatology  generated  over  Ft.  Devens  for  the 
month  of  August.  We  conclude  the  section  with  an  outline  of  the  operational  tests 
conducted  with  the  model  at  the  Turner  Drop  Zone  Tactical  Training  Area  (Turner 
DZ).  A  summary  of  each  individual  case  for  the  Turner  DZ  tests  is  presented  in 
Appendix  A.  We  recommend  that  operational  forecasters  read  Appendix  A  to  get 
acquainted  with  the  methods  we  used  to  determine  appropriate  input  conditions 
for  running  the  model  in  real  time  on  a  day-to-day  basis. 

5.1  Detachment  12  Mission  and  Forecaster  Training 

Detachment  12  is  essentially  a  two-faceted  operation.  The  base  weather 
station,  located  at  Moore  Army  Airfield  (Moore  AAF),  is  primarily  responsible 
for  preparing  the  TAF  and  briefing  aircrews  who  transit  the  area.  The  Special 
Operations  Weather  Team  (SOWT)  supports  Det  12's  primary  customer,  the  Army 
10th  Special  Forces  Group  (10  Group).  The  SOWT  deploys  with  the  10th  Group  on 
exercises  and  missions,  and  supports  paradrop  training  that  takes  place  at 
Turner  DZ,  located  11  km  southwest  of  Moore  AAF.  The  SOWT  consists  of  four 
jump-qualified  AWS  personnel;  three  forecasters  and  one  observer/forecaster 
apprentice. 

As  one  might  expect,  the  winds  measured  at  the  base  weather  station  (Moore 
AAF)  do  not  always  agree  with  the  winds  measured  at  Turner  DZ,  because  of  the 
11 -km  separation.  Another  problem  with  the  wind  measurements  at  Moore  AAF 
is  that,  due  to  terrain  effects,  windspeeds  for  directions  of  west -northwest  through 
north-northeast  are  consistently  underestimated  by  the  equipment  by  2  to  4  kt. 
Forecasters  typically  add  5  kt  to  the  TAF  wind  forecast  to  account  for  this  problem. 
In  Section  5.  2,  we  describe  our  efforts  to  simulate  this  windflow  using  the  model 
and  three-dimensional  depictions  of  the  terrain  surrounding  Moore  AAF  and 
Turner  DZ. 

The  weather  support  given  by  the  SOWT  to  the  10th  Group  consists  mainly  of 
weather  forecasts  and  observations  for  Turner  DZ,  and  climatologies  for  deploy¬ 
ment  locations  as  well  as  potential  "hot  spot"  areas  around  the  world.  The  10th 
Group's  paradrop  mission  is  perhaps  the  most  weather  sensitive  of  all  the  missions 
t he v  have.  They  are  particularly  vulnerable  to  winds;  surface  windspeeds  (including 


gusts)  =  13  kt  or  drop  altitude  [(900  to  1250  ft  AGL  for  personnel;  1700  ft  AGL  for 
cargo)]  windspeeds  ^  30  kt  essentially  curtail  operations  for  as  long  as  these  con¬ 
ditions  persist.  Wind  chill  effects  are  another  consideration  (chill  factors  ^  -25°F 
are  the  criteria  here).  Thus,  wind  forecasts  have  an  important  impact  on  10th 
Group's  ability  to  accomplish  its  mission.  In  particular,  the  windspeed  criteria 
presented  are  sensitive  values  that  are  not  easy  to  predict.  It  appears  that  the 
SOWT's  missions  are  well-suited  for  application  of  our  windflow  modeling  system. 

The  paradrop  operation's  sensitivity  to  the  wind  needs  to  be  further  discussed. 
Nearly  every  aspect  of  the  operation,  from  the  placement  of  marker  panels  over 
the  DZ,  to  the  well-being  of  the  troops  once  they  have  landed,  is  sensitive  to  the 

wind  direction  and  speed.  The  DZ  is  usuallv  set  up  over  an  open  field  area  approxi- 
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mately  1/2  km"  or  greater.  Marker  panels  are  placed  over  the  DZ  area  to  give 
the  aircraft  commander  a  landmark  he  can  use  to  guide  the  aircraft  into  the  DZ. 

The  markers  also  give  the  jumpmaster  (the  individual  responsible  for  releasing 
the  troops  from  the  aircraft)  an  idea  when  he  can  release  his  troops  to  begin  their 
.jumps.  The  timing  of  the  jumpmaster's  release  of  his  troops  is  affected  by  head¬ 
winds  or  tailwinds  over  the  DZ.  It  is  during  the  latter  part  of  the  jump  that  the 
winds  are  perhaps  most  crucial  when,  about  300-ft  AGL,  the  jumper  begins  to  turn 
his  bodv  into  the  wind.  then  a  light  wind  can  mean  the  difference  between  a  safe, 
soft  landing,  and  an  injurv  resulting  from  a  hard  landing  or  a  drift  into  the  treetops. 

Currently,  the  SOWT  has  several  means  for  estimating  winds  relative  to 
advising  go  or  no-go  decisions  for  a  jump.  The  first  way  is  through  the  TAF  wind 
forecast  that  provides  surface  wind  forecasts  up  to  24  h  before  a  jump. 

Another  wav  is  through  use  of  a  forecast  study  prepared  for  Ft.  Devens  that  uses 
a  nomogram  giving  the  surface  wind  speed  at  the  drop  zone  as  a  function  of  the 
pressure  gradient  and  the  wind  direction  derived  from  the  latest  surface  analysis 
(this  nomogram  is  shown  in  Figure  15).  The  third  way,  most  commonly  used,  is 
the  computation  of  a  layer-mean  wind  (surface  to  1 2  50  -  ft  AGL,  usually)  by  tracking 
a  10-g  balloon  launched  about  1-1.5  h  prior  to  t  ime -on -target  (TOT),  and  supple¬ 
mented  bv  surface  wind  observations  taken  at  the  DZ  by  the  DZ  Safety  Officer  using 
handheld  equipment.  The  laver-mean  wind  is  computed  through  knowledge  of  the 
balloon's  ascent  rate,  and  azimuth  and  elevat'on  angle  above  the  surface  as  mea¬ 
sured  with  a  scope.  None  of  these  methods  provide  wind  profile  data  at  different 
levels.  Although  the  model  wind  profiles  described  in  Section  4.3  only  extend  to 
300  ft  AGL,  there  are  three  points  we  believe  should  be  made  about  the  usefulness 
of  the  model  wind  profiles.  First,  recall  that  the  lowest  300  ft  are  the  most 
important  in  a  jump  (as  previously  explained).  Second,  the  drop  altitudes  in  war¬ 
time  deplovment  would  he  closer  to  500  ft  AGL.  Third,  future  parachute  equipment 


PRESSURE  GRADIENT  (MB/3deg  LAT) 


is  being  designed  for  safe  operational  drops  at  heights  as  low  as  300  ft  AGL 
57 

(Harclerode  ).  The  use  of  this  future  equipment  could  be  adequately  supported 
by  the  surface -layer  wind  profiles  produced  by  our  model. 
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WIND  DIRECTION 


Figure  15.  Nomogram  Describing  Wind  Forecast  Technique  for  Turner  DZ 
Operations.  A  three -degree  latitude  interval  is  oriented  over  surface 
analysis  until  strongest  pressure  gradient  over  the  station  is  identified. 
Windspeed  is  obtained  bv  reading  corresponding  pressure  gradient  value  and 
wind  direction  from  surface  analysis.  This  method  yields  the  interval  for 
maximum  surface  wind  gust  for  a  3  h  period  following  the  analysis  time 


Forecaster  training  on  the  windflow  modeling  system  consisted  of  six  sessions 
using  Dot  12's  Z-100  computer.  There  were  three  different  programs  on  which 
forecasters  were  trained  (see  also  Section  3).  The  first  one  was  the  windflow  model 
itself,  configured  to  run  interactively  with  the  user  (see  Figure  11  and  Section  4.2). 


57.  Harclerode,  P.  (l!*86)GQ's  8  metre  low-level  parachute.  Special  Force 
Journal  1:833-850. 


Forecasters  provided  input  data  either  from  current  observations,  the  TAF,  or 
most  frequently  the  TAF  updated  with  current  observations.  The  next  program 
was  a  horizontal  wind  plot,  that  automatically  read  in  the  output  files  from  the 
model.  The  plot  program  has  two  display  options:  winds  depicted  bv  arrows  whose 
length  is  proportional  to  speed,  or  the  standard  wind  plot  with  barbs.  An  example 
of  the  screen  display  for  the  first  DZ  case  (17  September  1986)  is  shown  in  Figure 
16.  A  terrain  contour  display,  also  shown,  was  provided  so  the  user  could  identify 

DAY,  MONTH,  HOUFM  7/SEP/1800 

MEAN  SPEED=5  KTS  MEAN  DIR=329  MEAN  EXP-.308 


STAB  PAR  =  20I  WIND  HEIGHT  IN  FEET=32.8 


Figure  16.  Model  Horizontal  Wind  Plot  From  Z-100  Sc  reen  for  17  September  1986 
at  1800Z.  Mean  wind  direction  and  speed  are  displayed  at  top,  as  is  mean  value 
of  model -generated  profile  exponent.  Stability  parameter  is  shown  at  bottom. 
Terrain  overlay  is  from  3X3  km  window  surrounding  Turner  I)Z 


geographic  reference  points  over  the  model  display  window.  The  third  program 
was  a  vertical  windspeed  profile  display  over  the  lowest  300-ft  AGL.  In  Section 
4.  3.  2  we  mentioned  that  the  model-generated  wind  profiles  tended  to  overestimate 
'he  winds  at  levels  above  the  IO-m  computation  height  (recall  Figure  13).  To  use 
the  profiles  for  real-time  support  of  Turner  DZ  operations,  we  developed  a  scheme 
based  on  calculating  the  wind  profiles  using  a  representative  height  of  31.6  m  (the 
geometric  mean  of  10  m  and  100  m).  This  method,  shown  in  Figure  17,  yields  more 
reasonable  values  of  the  wind  profile  exponent  than  those  based  on  the  model  run. 

To  illus'rate  this,  let  us  take  the  17  September  case  whose  wind  plot  is  shown  in 
Figure  16  and  derive  the  wind  exponent  for-  use  in  the  profile  program.  The  sta¬ 
le! .  •  v  class  from  Figure  16  is  about  2.01,  which  is  below  3.3;  therefore,  conditions 
arc  uns’able,  and  from  Figure  17,  the  exponent  is  0.  16.  The  model -generated 
exponent,  also  shown  in  Figure  16,  is  0.3  1,  which  is  almost  twice  the  derived  value. 
Plots  of  the  wind  profile  using  the  model -generated  exponent  and  the  derived  expo¬ 
nent  are  shown  in  Figures  18a  and  18b,  respectively.  A  comparison  of  the  two 
profiles  sh  iws  that  'he  modified  profile  is  flatter  than  that  generated  by  the  model, 
with  a  300-ft  windspeed  of  7  kt,  as  opposed  to  10  kt  from  the  model.  The  Mean 
Fffective  Wind  (MEW)  as  measured  for  this  case  was  310°/ 10  kt,  so  that  the  modi¬ 
fied  profile  is  slightly  more  consistent  with  the  Mean  Effective  Wind  than  is  the 
model -generated  profile.  For  all  of  the  13  jump  operations  for  which  we  provided 
model  support,  the  windspeed  profiles  based  on  Figure  19  were  more  consistent 
with  ground-  and  balloon -based  measurements  than  those  derived  from  the  model 
exponent.  We  should  note  here  that  none  of  our  13  cases  had  stable  conditions, 
and  the  profile's  reliability  is  still  suspect  under  these  conditions  (recall  Section 
4.  3.  2). 

All  of  the  personnel  whom  we  trained  agreed  that  the  windflow  modeling  system 
was  easv  to  learn  and  use,  and  that  the  system  was  a  valuable  addition  to  SOWT 
support  at  Turner  DZ. 

5.2  Model  Windflow  Climatology  Simulations  for  August 

In  order  to  demonstrate  the  uses  for  the  windflow  model  in  a  deployment 
scenario,  we  generated  windflow  climatology  over  the  Ft.  Devens  reservation, 
using  August  climate  data  from  the  Terminal  Forecast  Reference  Notebook  (T*FRN). 

1  he  TFRN  is  a  standard  AWS  publication  containing  station  climatology  and  local 
forecast  studies.  The  model  was  run  for  the  predominant  wind  directions  of  south¬ 
west  and  west  using  various  windspeeds  for  both  daytime  and  nighttime  conditions. 
These  wind  simulations  were  run  to  determine  the  presence  of  local  terrain  effects 
on  the  windflow,  particularly  over  Turner  DZ  and  Moore  AAF. 
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Figure  19.  Ft.  Devens  Model  Terrain  Elevations  in  m  MSL.  The  3  <  3  km 
window  used  lor  the  operational  model  test  is  in  the  west-rentral  portion 
of  the  model  domain.  Turner  D Z  is  highlighted  by  the  dashed  border  within 
the  window.  The  area  surrounding  Moore  AAF  is  shown  bv  solid  outline  of 
northeast  corner  of  domain 


The  model  domain  for  the  simulations  is  shown  in  Figure  19.  The  domain 
size  was  8.  8  km  west  to  east  by  13.  8  km  south  to  north.  With  a  200-m  horizontal 
spacing  between  grid  points,  this  translates  into  an  array  44  X  69  (3036  grid  points). 
Since  we  had  no  vegetation  height  data  for  this  area,  we  assumed  a  constant  sur¬ 
face  roughness  of  50  cm,  corresponding  to  gentle  topography  with  some  forest 
cover.  Turner  DZ  (shown  by  the  box  in  Figure  20)  is  located  in  a  large  open  area 
(900  X  600  m)  of  sloping  terrain.  The  area  is  well-exposed  to  winds  with  directions 
from  southerly  to  northwesterly.  A  small  hill  (about  115  m  above  mean  sea  level 
MSL)  appears  in  the  northeast  corner.  The  terrain  slopes  to  the  southwest,  from 
about  115  to  120  m  MSL  to  around  90  m  MSL  in  the  southwest  corner.  The  average 
terrain  grade  is  about  1.7  percent.  Vegetation  over  the  I)Z  consists  of  tall  grasses 
in  a  sandy  soil.  The  DZ  is  surrounded  on  all  sides  bv  mixed  coniferous  and 
deciduous  trees  about  10  to  15  m  high,  Moore  AAF  (shown  in  Figure  21)  sits  on 
a  small  plateau  oriented  northwest  to  southeast,  about  1.7  km  long  and  0.8  km  wide. 
The  plateau  is  about  15  to  20  m  above  the  surrounding  terrain  that  consists  mainly 
of  forest,  several  marshes,  and  the  nearby  Nashua  River. 


Figure  20.  Enlargement  of  3  X  3  km  Window  Used  for 
Model  Tests  at  Turner  DZ.  The  DZ  is  highlighted  by 
the  solid  border  in  the  center 


Figure  2  1.  U.  S. 
Geological  Survey 
Map  at  1:25,000  Scale 
Showing  Terrain 
Elevations  (in  ft  MSL) 
Over  Area  Surrounding 
Moore  AAF.  The 
"G"  denotes  the 
location  of  GMQ-20 
wind -measuring 
equipment  at  the  north 
end  of  the  runway. 

The  "S"  denotes  the 
location  of  a  windsock 
some  300  m  south- 
southeast  of  ”G”. 

The  dashed  line  shows 
the  location  over 
which  a  cross  section 
(Figure  22)  was  taken 
to  TIustrate  the  flow 
separation/eddy  effect 


Recall  from  Section  5.  1  that  Moore  AAF  has  a  wind  equipment  problem  because 
of  the  terrain.  We  examined  the  terrain  around  Moore  AAF  using  horizontal  maps 
and  three-dimensional  terrain  diagrams  generated  from  our  terrain  elevation  data 
for  Ft.  Devens.  From  forecaster  interviews  and  the  TFRN,  we  determined  that 
the  w tnd -measurement  problem  was  being  caused  by  flow  separation  taking  place 
as  the  wind  encounters  the  sudden  terrain  slope  at  the  north  end  of  the  airfield. 

This  flow  separation  was  producing  small  turbulent  eddies  over  the  area  where  the 
wind  equipment  was  located  (see  terrain  cross  section  and  flow  diagram  in  Fig¬ 
ure  22).  However,  after  the  flow  rides  up  the  plateau's  side,  it  descends  about 
150  to  200  m  downwind  of  the  airfield  edge.  This  descent  has  been  confirmed  by 
hand-held  equipment  and  a  wind  sock  that  is  located  about  300  m  downwind  of  the 
wind  equipment.  Flow  separation  along  the  edges  of  escarpments  and  cliffs  has 

been  studied  by  several  researchers,  such  as  Scorer"^’ ^  ancj  Jones.  A  con- 

6 1 

cise  summary  of  flow  separation  and  wake  effects  appears  in  Orgill.  Observa¬ 
tions  of  flow  separation  show  that  the  zones  of  turbulent  eddies  at  the  edge  of  a 
cliff  c  an  be  2  to  10 -m  deep  (well  above  the  level  of  a  standard  wind  measurement 
system),  and  can  extend  about  50  to  200-m  downwind.  Thus,  the  observations  at 
Moore  AAF  are  confirmed  by  other  observational  studies.  The  horizontal  scale 
of  these  eddies  is  too  small  to  be  accounted  for  by  the  windflow  model.  Numerical 
models  can  only  "recognize"  features  having  a  wavelength  equal  to  or  greater  than 
four  grid  spaces.  Our  model,  with  a  200-m  spacing,  can  only  resolve  features 
with  wavelengths  800  m  or  greater.  Another  way  of  phrasing  this  would  be  to  say 
that  if  we  wanted  to  attempt  to  model  this  flow  separation,  we  would  need  a  grid 
spacing  of  about  50  m  (this,  of  course,  does  not  guarantee  that  we  could  simulate 
this  feature  simply  by  reducing  the  grid  spacing). 

58.  Scorer,  R.S.  (1955)  Theory  of  airflow  over  mountains:  IV -Separation  of  flow 

from  the  surface,  Quart.  J.  Roy.  Meteorol,  Soc.  8 1  ;340-350. 

59.  Scorer,  R.S.  (1972)  Clouds  of  the  World,  Stackpole  Books,  Harrisburg,  PA. 

60.  Jones,  O.  K.  (1970)  The  flow  of  a  stratified  fluid  over  a  vertical  step,  Tellus 

22(No.  5);48  1  -492. 

61.  Orgill,  M.  M.  (1981)  Atmospheric  Studies  in  Complex  Terrain;  A  Planning 

Guide  for  Future  Studies,  PNL-3656,  ASCOT/80/4,  pacific  Northwest 


1-' i tin t-f  22.  Illustration  of  Flow  Separation/Kddv  Rffeet  Occurring  With  a  Wind 
Direction  nt  240  ,  The  cross  section  shows  the  terrain  elevations  in  m  with 
'he  s  t  ream  I  mes  and  arrows  describing  the  flow  up  to  about  .20  m  AG  L.  The 
GMQ-20  equipment  is  depicted  bv  the  windvane;  the  windsock  is  shown  about 
200  m  downwind 

Although  the  model  cannot  account  for  differences  in  the  wind  measurements 
along  'he  runwav  of  Moore  AAF,  it  can  give  us  information  about  slope  wind  effects 
on  the  flow.  We  accomplished  this  bv  comparing  the  mean  wind  directions  and 
speeds  for  Moore  AA  F  and  Turner  1)7.  for  each  set  of  climatological  simulations. 

As  these  mean  winds  get  closer  together,  the  terrain  effects  can  be  said  to  have 
decreased  over  the  Ft.  Devons  region.  We  compared  Moore  AAF  to  Turner  DZ 
to  see  if  :inv  svstematic  wind  differences  between  the  two  locations  could  be  found. 
This  information  would  be  useful  in  that  a  relocation  of  the  wind  equipment  at 
Moore  AAF  would  not  necessarily  eliminate  discrepancies  in  wind  measurements 
between  the  airfield  and  the  DZ.  The  results  of  our  wind  microclimate  simulta- 
tions  appear  in  Table  6.  We  are  not  surprised  to  see  that  slope  wind  effects  are 
mlv  important  for  windspeeds  below  6  kt,  and  are  essentially  negligible  above  6  kt. 
Helow  6  kt,  the  winds  at  Turner  DZ  are  about  1-2  kt  stronger  than  at  Moore  AAF’ 
luring  the  dav,  but  at  night  the  speeds  are  virtually  the  same.  Wind  directions 
at  Turner  DZ  also  'ended  to  be  about  10  to  35’  more  southerly  than  at  Moore  AAF’ 
during  the  dav.  At  night,  winds  at  Turner  were  5  to  30’  more  northwesterly  than 
at  the  airfield.  Keep  in  mind  that  these  simulation  results  are  only  applicable  for 
sou t hw'i's tor! v  and  westerly  winds  at  Moore  AAF"  and  for  essentially  clear  to 
scattered  conditions.  For  cloudv,  neutral  stability  conditions,  we  would  expect  the 
differences  between  the  two  locations  to  be  smaller.  From  Table  6  we  can  conclude 


Table  6.  Summary  of  August  Wind  Climatology  Simulations  for  Ft.  Devens 

(August  15) 


Daytime  (1400  LST ) 

Input  Conditions:  Temperature  =  80°  F.  Cloud  Cover  =  1/8 


Input  Wind 

Di  r  'Speed 

No.  Step 
to  Solution 

Model  Mean  Dir/Speed  (Max  Speed) 

Moore  AAF  Turner  l)Z 

270’1  kt 

26 

279’/  1 

kt  (2.  4) 

2 4 5 ° / 2  kt  (2. 7) 

27  0  V'4  kt 

23 

273  ’/3 

kt  (4.  0) 

265° / 5  kt  (5.0) 

27  0  V'  6  kt 

5 

27  2  °  /  6 

kt  (6.  1) 

269 ’/ 6  kt  (6. 3) 

27  0’  '8  kt 

2 

27  l’/8 

kt  (8.  1) 

272 ’ /8  kt  (8.  1) 

27  0’  /  16  kt 

2 

27  1  ’  /  1 6  kt  (16.  1) 

271’/  16  kt  (16.3) 

225  ’/ 4  kt 

34 

229  ’  / 3 

kt  (3. 4) 

216°/4  kt  (4. 2) 

225V'  6  kt 

15 

2  2  8  ’  /  5 

kt  (5.4) 

2 2 6 0  /  5  kt  (5.6) 

225'/ 8  kt 

2 

226°  / 8 

kt  (8.  0) 

226’/8  kt  (8.  1) 

225’/  16  kt 

2 

226 0 /  16  kt  (16. 2) 

226 °  / 16  kt  (16. 2) 

Nighttime  (2100  LST) 

Input  Conditions:  Temperature  =  70'  F,  Cloud  Cover  =  1/8 


270’/  1 

kt 

18 

263’/  1 

kt  (2.  2) 

295°  /  1 

kt  (2.  2) 

27  0  ’  /  4 

kt 

10 

268 '/ 4 

kt  (5.  2) 

279’  /4 

kt  (5.  2) 

270 '/ 6 

kt 

o 

27  l’/6 

kt  (6.  1) 

27  2’  16 

kt  (6.  1) 

270'  / 8 

kt 

2 

27  1  ’  /  8 

kt  (8.  2) 

270 ’/  8 

kt  (8.  2) 

270’/ 16  kt 

2 

27  1’/ 16  kt  (16.  1) 

271/16  kt  ( 16.  3) 

225’/ 4 

kt 

20 

220’ / 4 

kt  (5.  2) 

227  ’/  3 

kt  (5.  2) 

225  ’/  6 

kt 

7 

228 ’ / 6 

kt  (6.  0) 

2  2  9  ’  /  6 

kt  (6.  0) 

225’  /8 

kt 

2 

226’/8 

kt  (8.  0) 

226  '  /8 

kt  (8.  0) 

225’  /  16 

kt 

2 

226’/ 16  kt  ( 16. 2) 

226’/ 16  kt  (16. 2) 

that,  according  to  the  model,  terrain  slope  effects 

on  the  wdnds 

can  add  1 

-2  kt  to 

the  wind  ; 

speed  differences 

between 

Turner  I)Z  and  Moore  AAF. 

Adding  the  differ- 

enc  es  from  wind  sensor  plai  ement  on  the  runway  to  this,  we  fine]  that  a  reported 
light  wind  a!  Moore  AAF  mav  cause  problems  lor  jumpers  at  the  l)Z  as  illustrated 
in  'he  following  example. 


Daytime  (Unstable)  Conditions  with  Scattered  Clouds 


Moore  AAF  observation:  27(W5  kt 

Turner  DZ  winds: 

(1)  Add  slope  wind  effect  (-10  to  -35  =  /+l  to  2  kt):  235  to  260:>/6-7  kt 

(2)  Add  wind  instrument  error  (+2  to  4  kt)  :  235  to  260°/8-ll  kt 

(3)  Compute  wind  profile  exponent  (Figure  17)  :  Exponent  =  0.  16 

(4)  Profile  winds  at  300  ft:  11  to  16  ft  (critical  windspeeds  could  be 
encountered). 

The  preceding  example,  though  hypothetical,  illustrates  the  usefulness  of  the  model 
windflow  climatology  together  with  other  known  information  to  make  a  determina¬ 
tion  of  safety  conditions  at  the  DZ.  Situations  similar  to  the  above  example  have 
actually  been  observed  by  SOWT  members  at  the  DZ  on  several  occasions.  A  few 
of  our  cases  in  the  next  section  also  illustrate  these  effects,  using  both  the  wind- 
flow  model  output  and  vertical  wind  profiles  in  support  of  jump  operations  at 
Turner  DZ. 

5.3  Operational  Testing  for  Paradrop  Operations 
at  Turner  Drop  Zone 

This  section  summarizes  the  operational  tests  over  Turner  DZ  using  the  model 
and  its  peripheral  programs.  We  followed  the  guidelines  for  model  input  that  were 
outlined  in  Sections  4.2  and  4.3.  The  model  domain  used  for  this  real-time  test 
was  a  3  X  3  km  window  (highlighted  by  solid  lines  in  Figure  19;  window  itself  shown 
in  Figure  20)  surrounding  Turner  DZ.  The  window  was  a  15  X  15  array  (225  grid 
points),  and  was  constrained  by  the  limitations  of  Det  12's  Z-100  computer  without 
the  "math  chip.  "  As  a  rule,  the  model  took  12  to  20  min  to  produce  a  simulation. 
Use  of  TAF  information  for  input  enabled  the  model  to  be  run  as  a  forecast  model, 
giving  SOWT  forecasters  as  much  as  24  h  lead  time  in  examining  horizontal  winds 
and  vertical  profiles  for  a  jump  operation. 

Operational  testing  was  usually  conducted  in  the  following  manner.  A  pre¬ 
liminary  ("first  look")  model  run  was  made  as  much  as  24  h  prior  to  a  jump  using 
input  from  the  TAF.  On  the  morning  before  the  scheduled  jump  (most  jumps  were 
in  the  afternoon,  but  we  had  two  night  jumps),  we  reran  the  model  using  the  latest 
TAF  updated  with  current  observations.  We  sometimes  ran  the  model  more  than 
once  using  the  most  current  observations.  Approximately  3  to  6  h  prior  to  TOT  we 
chose  the  model  run  we  believed  was  most  accurate  and  used  it  for  supporting  that 
day's  operation.  This  method  allowed  us  to  begin  looking  at  forecast  conditions 
over  the  DZ  the  day  before  a  jump,  and  enabled  us  to  "metwatch"  DZ  conditions 
right  up  until  the  time  of  balloon  launch  1-1.5  h  before  TOT. 


The  summary  of  our  operational  testing  is  shown  in  Table  7.  There  are 
several  interesting  results  from  the  evaluation.  First,  notice  that  the  model  wind 
profile  indicated  300  ft  winds  in  excess  of  the  critical  surface  value  on  five  occa¬ 
sions.  On  four  of  these  occasions  wind  gusts  greater  than  13  kt  were  reported  at 
•  he  DZ  and  on  a  fifth  occasion  (28  October  1500Z),  gusts  in  excess  of  critical  speed 
appeared  about  1.  5  h  after  verification.  The  model  wind  profiles  showed  these 
st  rone  winds  aloft  despite  the  fact  that  neither  the  input  TAF  nor  the  model  surface 
wind  plot  indicated  windspeeds  d  13  kt.  We  believe  that  the  model  wind  profiles 
can  he  used  to  predict  gustv  surface  winds  over  the  DZ  because  the  winds  from  aloft 
can  be  mixed  down  to  the  surface  through  turbulent  eddy  processes  (recall  Section 
4.  lh  Another  interesting  result  of  our  evaluation  is  that  the  model  correctly  diag¬ 
nose  i  -he  stabi lit v  as  neutral  instead  of  stable,  for  the  two  nighttime  cases  (2  and  3) 
•hat  had  st  ratocumulus  cloud  decks.  This  correct  stability  diagnosis  led  to  model 
wind  profiles  that  agreed  well  with  the  measured  mean  effective  winds.  Notice 
trim  Table  7  ’hat  there  were  also  several  occasions  (cases  4,  7,  9)  that  behaved 
similarlv  to  our  example  presented  at  the  end  of  the  last  section.  The  most  dra¬ 
matic  case  was  No.  4,  where  a  10  kt  difference  was  observed  between  Moore  AAF 
and  Turner  DZ.  The  wind  direction  at  the  DZ  was  50°  less  than  at  the  airfield,  and 
surface  gusts  in  excess  of  the  critical  value  were  observed  at  the  DZ.  In  all  three 
cases,  the  wind  trajectories  were  from  the  west  to  northwest  and  were  consistent 
with  the  flow  separation  theory  discussed  in  the  previous  section  and  described  in 
die  schematic  diagram  of  Figure  22.  This  last  result  illustrates  the  benefits  from 
producing  model  wind  climatologies  and  combining  the  results  with  other  known 
information  about  the  station. 

Although  the  results  documented  in  Table  7  look  impressive,  we  should  make 
several  cautionary  points  about  our  study.  First,  this  is  a  small  sample  size  (13 
■  a  ,i's)  that  consisted  mainly  of  neutral  cases.  For  stable  cases,  we  advise  caution 
in  'lie  interpretation  of  the  model  profiles.  We  would  expect  some  model  "false 
alarms"  to  appear  in  a  larger  data  sample.  We  must  also  remember  that  the  model 
output  is  only  as  good  as  the  TAF  or  updated  information  that  is  put  into  the  model. 
Notice  from  Table  7  that  the  greatest  error  between  input  and  observed  DZ  wind- 
speeds  was  only  about  5  kt  (cases  2  and  12).  This  was  due  to  the  quality  TAFs  we 
received  and  the  diligent  updates  that  were  made  before  running  the  model. 
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Appendix  A 

Description  of  Tumor  Drop  Zoeo  Coco  StoAoo 


This  appendix  is  structured  in  the  following  manner; 

(1)  Summary  of  1200Z  synoptic  situation  for  each  case,  to  include 
relevant  surface  features  (for  example,  fronts,  highs,  lows, 

and  so  on),  predominant  surface  and  upper-level  (500-mb)  flows, 
with  information  about  precipitation  during  the  previous  24  h. 

(2)  Discussion  of  input  to  model  (cloud  cover /type/ level,  visibility, 
temperature,  winds),  plus  any  relevant  information  concerning 
multiple  model  runs,  choice  of  runs,  and  so  on. 

(3)  Discussion  of  conditions  at  Turner  DZ  for  each  operational  case. 

For  Cases  2  and  3,  the  synoptic  discussion  covers  the  period  from  the  pre¬ 
ceding  1200Z  analysis  to  the  following  1200Z  analysis. 

Case  1  —  17  September  1800Z 

Synoptic  Situation  -  High-pressure  ridge  from  northern  Quebec  to  mid-Atlantic 
states.  Light  northwesterly  flow  over  New  England.  Dry  northwest  flow  over  the 
northeast  at  500  mb.  Light  precipitation  recorded  over  eastern  half  of  Massachu¬ 
setts  from  previous  24  h. 

Input  Conditions  -  Used  11Z  TAF  to  initialize  model.  Forecast  from  TAF: 

2/8  Sc  050,  7  +  mi  vsby,  65°F,  330° /5  kt. 
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DZ  Conditions  -  Paradrop  may  by  C-130  aircraft  carrying  about  10  troops. 
Flight  level  was  about  800  to  1000-ft  AGL,  aircraft  approach  from  289°  (meteorol¬ 
ogical  direction).  Jumpers  were  taken  by  wind  towards  the  southeast,  and  one 
lumper  landed  close  to  trees  at  eastern  edge  of  DZ. 

Case  2—17  October  0100Z 

Synoptic  Situation  (Using  16  October  and  17  October  charts)  -  Weakening  cold 
front  move  1  from  northern  plains  and  Great  Lakes  to  Ohio  River  Valley.  Associa¬ 
te  1  low  and  cold  front  moved  from  north  of  Lake  Huron  to  western  Pennsylvania 
in  i  New  York  Citv  area.  High  pressure  located  over  Ontario  -  Manitoba  border 
moved  east  to  Lake  Superior  and  western  Quebec.  Light  and  variable  flow  over 
New  Fngl.m  1  became  northerly  after  frontal  passage.  Trough  at  500-mb  centered 
i bout  Hf>  VV  shifte  i  to  about  ROW  as  several  short  waves  moved  through  flow.  No 
!>r>  c l  m t a •  ion  luring  period  over  Massachusetts. 

Input  Coalitions  -  Used  17Z  16  October  TAF  to  initialize  model.  Forecast 
conditions;  a '8  Sc  035  F/8  Ac  100  8/8  Cs  250,  7  f  mi  vsby,  46’  F,  290'/5  kt. 

I  >/.  Conditions  -  Balloon  ascent  showed  direction  about  295’  over  lowest  several 
bun  Ire  I  feet  before  beginning  to  drift  towards  the  south.  Paradrop  from  two  C-130s 
living  between  1000  ft  and  1100  ft  AFL.  About  four  runs  were  made  by  each  air- 
crib,  Iropning  about  eight  to  nine  troops  per  run.  Several  different  approach 
lire,  riotis  (from  250  ’  to  205  )  produced  cross-wind  components  to  jumpers  causing 
both  'he  marker  panels  and  aircraft  paths  to  be  adjusted  after  the  first  run  (when 
several  ■  mops  drifted  too  close  to  eastern  edge  of  DZ).  Drops  were  made  between 
about  004 0  an  I  0130Z.  Most  jumpers  experienced  a  southward  drift  due  to  pre¬ 
vailing  win  Is  (see  Table  8). 

Case  3  —  Iff  October  0630/ 

Svnop'ir  Situation  (Using  18  October  12007.  charts)  -  High  pressure  centered 
over  Lake  Huron.  Northerly  to  northwesterly  flow  over  New  England.  Sharp 
500-mb  trough  over  New  England.  Light  pre<  ipitat ion  fell  over  much  of  northeast 
during  previous  24  hr. 

Input  (  ,,n  iitmns  -  Used  17/ Oct  TAF  for  input.  Forecast  conditions:  7/8  Sc 
035  4  Ac  080,  7  »  mi  vsbv.  38  F,  010  / 5  kt. 

It/  Conditions  -  Balloon  launch  showed  an  almost  vertical  rise  due  to  light 
winds.  Wind  hreeijon  was  uniform  from  north -northeast  most  of  way,  turning  to 
northeast  near  top  of  Inver.  Paradrop  from  C-130  flying  at  standard  drop  altitude. 
Am  raft  made  live  runs  across  l)Z  before  dropping  troops  on  sixth  run.  About  six 
drops  were  m  i d» ■  from  065D  to  0730/..  Thirty-eight  National  Guardsmen  made  the 
lump,  and  despite  the  light  winds,  several  who  rode  with  the  wind  (recall  Figure  16) 
came  down  unite  rapidly. 


Case  4  —  24  October  1800Z 


Synoptic  Situation  -  Series  of  high  pressure  cells  extended  from  British 
Columbia  through  Ontario.  Weak  stationary  front  from  southern  Iowa  eastward 
to  Pennsylvania,  becoming  cold  front  off  Delaware  Peninsula.  Light  northwesterly 
flow  over  New  England.  Strong  cyclonic  flow  at  500  mb  around  trough  located  off 
New  England  coast.  Precipitation  had  occurred  over  Massachusetts  during  period. 

Input  Conditions  -  Used  11Z  TAF  to  initialize  model.  Forecast  conditions: 

1/8  Ac  100  2/8  Ci  250,  7  +  mi  vsby,  57  3  F  350  /5  kt,  wind  350=  /  10-  15  kt  14-00Z. 
Model  runs  made  first  for  350’/5  kt.  then  for  350’/ 10  kt  (from  TAF  remark).  We 
decided  to  use  second  model  run  for  verification  on  the  basis  of  5-6  kt  winds 
reported  at  Moore  AAF  that  morning. 

DZ  Conditions  -  Balloon  showed  strong  winds  between  100'  and  500'  AGL, 
then  decreasing  above  500'  to  produce  MEW  of  320’/ 12  kt.  Personnel  at  DZ 
(SOWT  and  DZSO)  have  seen  this  type  of  speed  profile  before;  however,  no  problems 
occurred  during  course  of  the  jump  from  C-130  of  approximately  60  troops  from 
900  to  1000  ft  AGL.  It  is  possible  that  strong  low-level  winds  recorded  by  balloon 
could  have  been  due  to  transient  turbulent  eddies,  as  evidenced  by  the  surface  gusti- 
ness  of  the  wind  this  day.  Jumps  were  made  from  1815  to  1900Z,  and  again  from 
2000  to  2100Z.  At  2000Z,  Moore  AAF  recorded  winds  of  360’ /6  kt  with  gusts  to  13  kt. 

Case  5  —  26  October  1500Z 

Synoptic  Situation  -  Cold  front  over  northern  New  England,  developing  warm 
front  over  mid-Atlantic  states,  with  light  rain  moving  into  southern  New  England 
from  southwest.  Very  light  winds  over  the  region.  Ridge  at  500  mb  over  New  York 
and  New  Jersey  with  westerly  flow  over  New  England. 

Input  Conditions  -  Since  no  TAF  was  available  on  the  weekend,  we  used  a 
combination  of  current  regional  observations  to  initialize  the  model.  Aware  of  the 
wind  problems  at  the  airfield  (they  were  open  for  observations),  we  surveyed  nearbv 
stations  and  found  them  all  to  have  light  winds,  with  ORB  (Worcester,  Massachusetts) 
reporting  a  direction  of  350  .  Input  conditions  were  as  follows:  8/8  As  100, 

6  mi  vsby.  light  rain  and  fog,  46  F,  350  /  I  kt. 

DZ  Conditions  -  Balloon  launch  indicated  verv  little  wind  until  2000  to  3000  ft 
AGL.  Paradrop  of  about  40  troops  from  C-130  took  about  45  minutes.  One  injury 
occurred,  but  was  not  weathe  r  -  re  lated. 

Cases  6  and  7  -  28  Oc  tober  1300Z.  1500Z 

Synoptic  Situation  -  Cold  front  was  located  about  500  km  offshore.  Weak  high 
pressure  ridge  from  Quebec  to  Ohio  River  Valley.  Light  and  variable  flow  over 
New  England.  Trough  at  500  mb  from  New  York  to  North  Carolina,  and  a  jet  streak 
moving  over  New  England.  Light  precipitation  had  fallen  over  northeastern  third  of 
1 '.  S.  during  pe riod. 
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Input  Conditions  -  Used  18Z  27  October  TAF  for  input.  Runs  were  made  pre¬ 
vious  day.  Forecast  conditions  prior  to  13Z;  8/8  St  012,  1/2  mi  vsby  in  fog, 

519F,  300’/5  kt.  Conditions  after  13Z:  4/8  St  012  5/8  Sc  030,  3  mi  vsby  in  fog, 
51”F,  3309  /8  kt.  Model  runs  were  made  using  conditions  prior  to,  and  after  13Z, 
respectively.  The  first  run  was  used  for  13Z  verification,  while  the  second  run 
was  made  for  15Z  verification. 

DZ  Conditions.  First  balloon  launch  at  1300Z  showed  light  winds  until 
200  to  300  ft,  when  it  began  to  accelerate,  indicating  higher  windspeeds  from  300  ft 
to  flight  altitude.  Paradrop  from  Army  Rlackhawk  helicopter  was  delaved  from  137. 
to  15Z  because  of  fog  at  Moore  AAF.  The  balloon  launched  at  157.  showed  a  more 
uniform  wind  distribution  with  height.  Paradrop  of  about  42  troops  lasted  from 
1520  to  1645Z.  No  problems  during  operation  until  last  jump  when  one  jumper 
landed  in  trees  because  he  did  not  turn  into  the  wind  soon  enough.  A  second  jumper 
landed  on  the  dirt  road  close  to  I)Z  entrance.  By  the  last  two  tumps,  the  surface 
winds  had  increased  to  12  kt  with  stronger  gusts.  We  suggest  that  strong  vertical 
mixing  was  beginning  to  bring  down  stronger  winds  from  aloft  at  this  time.  This 
observation  is  supported  "bv  the  13Z  and  157.  balloon  tracks,  that  changed  from 
non-uniform  to  uniform  wdnd  conditions,  and  bv  the  appearance  of  a  broken  Sc  dec  k 
at  3000  ft  at  17  Z. 

Case  8-29  October  1900Z 

Synoptic  Situation  -  High  pressure  off  V  irginia  coast.  Cold  front  approaching 
from  Great  Lakes.  Southwesterly  flow  in  advance  of  front  over  New  England. 
Southwest  flow  at  500  mb  over  New  England  through  flat  ridge  just  off  coast.  Spotty 
precipitation  over  Massachusetts  during  previous  period. 

Input  Conditions  -  Two  model  runs  were  made,  using  the  187.  28  October  and 
12Z.  29  October  TAFs,  respectively.  Forecast  conditions  from  first  TAF:  4/8  Sc. 

7  *  mi  vsby,  68  '  F,  240’  /  10  kt.  Forecast  from  second  TAF:  2/8  Sc  040  2/8  Ci  250, 
7  *  mi  vsby,  68 '  F,  270  \/6  kt.  We  chose  to  use  the  first  model  run  on  the  basis  of 
the  winds  observed  over  the  region  that  morning  (around  16Z),  which  were  south¬ 
westerly  with  speeds  of  7  to  15  kt. 

1)7.  Conditions  -  Balloon  track  showed  direction  about  2  10  below  300  ft,  shifting 
to  250'  above  300  ft  with  a  steady  increase  in  speed.  We  computed  MEWs  for  the 
following  levels:  0  to  500  ft,  250  /  1  1  kt;  0  to  1000  ft,  250"  /  13  kt;  0  to  1250  ft, 

250’/ 14-15  kt.  Gusty  surface  winds  were  observed  over  the  1)7.  at  this  time,  .lump 
scheduled  from  C-7  aircraft  was  cancelled  due  to  mechanical  problems.  We  be¬ 
lieve  it  mav  have  been  just  as  well  considering  the  strong  wind  conditions  at  the  DZ. 


Cases  9  and  10  —  30  October  1430Z,  1630Z 

Svnoptic  Situation  -  Weakening  cold  front  was  just  passing  through  the  Ft.  Devens 
area  at  1200Z.  Northwest  flow  was  behind  the  front.  Strong  west  to  southwest  flow 
at  500  mb  accompanied  this  system.  Precipitation  was  confined  to  New  York  and 
northern  New  England  with  this  system. 

Input  Conditions  -  We  made  two  model  runs  using  the  12Z  TAF  for  this  case. 

The  first  run  used  pre-frontal  conditions;  the  second  run  used  post-frontal  condi¬ 
tions  (at  17Z).  Pre-frontal  forecast:  2/8  Sc  035,  7  +  mi  vsby,  59’F,  240  to  270c>/4  kt. 
Post-frontal  forecast:  2/8  Sc  035,  7  +  mi  vsby,  60°F,  290°/ 12  kt.  It  was  obvious 
from  the  surface  analysis  that  the  frontal  passage  had  already  taken  place  by  12-13Z. 
We  therefore  used  the  post-frontal  run  (time  of  17Z)  for  verification.  The  model 
horizontal  wind  plot  for  this  case  showed  isolated  areas  of  15  kt  winds  over  the  DZ. 

DZ  Conditions  -  Balloon  launch  at  34  15Z  showed  steady  low-level  winds  be¬ 
coming  lighter  above  1000  ft  AGL.  Gusty  surface  winds  (peak  gust  of  19  kt  at  1530Z) 
hindered  Blackhawk  jump  operations  for  about  45  troops  from  1520  to  1630Z,  even 
causing  several  passes  to  be  held  until  the  gusts  subsided.  A  second  balloon  was 
launched  at  1640Z  by  Rhode  Island  Air  National  Guardsmen  for  C-130  cargo/equip¬ 
ment  drop.  Mean  effective  winds  were  calculated  at  different  levels  for  this  launch: 

0  to  600  ft,  320’  /  9  kt;  0  to  800  ft,  320’ /9  kt,  0  to  1100  ft,  320°/ 12  kt.  Winds  at 
Moore  AAF  were  330° /7  kt  with  gusts  to  12  kt  at  16Z. 

Case  11—3  November  1800Z 

Synoptic  Situation  -  High  pressure  centered  over  northern  New  England.  A 
ridge  of  high  pressure  extends  southwestward  into  the  southern  Great  Plains.  An 
approaching  frontal  system  >s  moving  into  the  Great  Lakes.  Light  and  variable 
surface  flow  exists  over  New  England.  A  strong  westerly  jet  at  500  mb  extends  from 
the  Great  Lakes  eastward  into  the  Canadian  Maritime  Provinces.  Light  precipita¬ 
tion  has  fallen  over  all  of  Massachusetts  during  the  previous  24  h. 

Input  Conditions  -  Used  12Z  TAF  to  initialize  model.  Forecast  conditions: 

1/8  Sc  050,  7  +  mi  vsby,  51’ F,  180’/7  kt  after  16Z. 

DZ  Conditions  -  Balloon  showed  fairly  moderate  winds  through  the  lowest  layers, 
becoming  lighter  above  700  to  750  ft.  Paradrop  of  eight  troops  from  C-130  took 
place  around  1835Z,  with  no  problems  encountered. 
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Case  12—5  November  1800Z 

Synoptic  Situation  -  Strong  high  pressure  over  the  Canadian  border  area  of 
New  England.  The  ridge  extends  southward  into  the  mid-Atlantic  states,  where  a 
stationary  front  lies  east-west  from  Cape  Hatteras,  North  Carolina  to  a  low  in 
northwestern  Mississippi.  Overrunning  precipitation  associated  with  the  stationary 
front  extends  from  the  Mississippi  River  Valley  eastward  into  central  New  Jersey. 
Surface  flow  over  New  England  is  northerly,  and  the  500-mb  flow  is  from  the  west- 
southwest  ahead  of  a  short  -wave  trough  in  Missouri.  Light  precipitation  has  fallen 
over  Massachusetts  during  the  previous  24  h. 

Input  Conditions  -  Used  12Z  TAF  to  initiate  model.  Forecast  conditions 
from  TAF;  5/8  Ac  080  7  Cs  250,  7  +  mi  vsby,  39°  F,  040<’/7  kts;  17-18Z:  5/8  Ns  030, 
5  mi  vsbv  in  light  fog  and  rainshowers,  35^F,  090°/7  kt;  19-20Z;  8/8  Ns  010,  3  mi 
vsbv  in  light  intermittent  rain,  35’ F,  100^/10  kt.  Information  from  the  TAF  for 
18Z  and  later  was  used  in  determining  input  conditions.  We  made  two  model  runs: 
one  using  the  pre-17Z  conditions,  and  the  second  using  the  17-18Z  conditions. 
Indications  in  the  late  morning  (15Z)  were  that  the  pre-17Z  conditions  would  remain 
over  the  area  at  jump  time,  so  this  forecast  was  used  for  verification. 

DZ  Conditions  -  Balloon  track  showed  light  northeast  winds  below  200  to  250  ft, 
changing  to  light  easterly  winds  above.  Paradrop  from  C-130  of  about  40  troops 
occurred  between  1820Z  and  1850Z.  Wind  drift  of  troops  was  slight,  even  for  those 
riding  with  the  winds. 

Case  13—7  November  1900Z 

Synoptic  Situation  -  Strong  high  pressure  is  centered  over  the  Connecticut  shore. 
The  ridge  extends  from  Nova  Scotia  into  the  southern  Appalachian  mountains.  Sur¬ 
face  flow  is  light  and  variable  over  New  England,  and  west  to  southwesterly  flow 
at  500  mb  over  the  northeast  exists  in  an  area  of  weak  ridging  from  New  England 
southward  to  the  mid-Atlantic  states.  Some  substantial  precipitation  amounts  (in 
excess  of  0.  25  in.  )  were  recorded  over  Massachusetts  during  the  previous  24  h, 
and  snow  fell  through  northeastern  New  York  and  most  of  Vermont,  New  Hampshire 
and  Maine. 

Input  Conditions  -  Used  12Z  TAF  to  initialize  model.  Forecast  conditions: 

2  '8  St  010  7/8  St  023,  5  mi  vsby  in  light  fog,  50°  F,  variable  winds  at  5  kt;  17 -18Z: 

4/8  Sc  045,  7  ♦  mi  vsbv,  57’F,  150'/7  kt.  Model  runs  were  made  from  both  fore¬ 
casts,  using  an  18Z  time.  In  late  morning  (around  15Z),  we  decided  to  use  pre-17Z 
forecast  conditions  for  verification.  We  used  a  220°  wind  direction  for  the  first 
run  based  on  observations  at  Worcester. 
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DZ  Conditions  -  Balloon  track  showed  southwest  winds  through  the  lowest 
350  ft,  becoming  west -southw-est  above  this  level.  Wind  speeds  remained  fairly 
uniform  throughout  the  layer.  Paradrop  from  C-130  of  approximately  30  para¬ 
troopers  took  place  between  1900  and  1945Z  (original  time  was  supposed  to  be 
1800Z).  Winds  were  not  a  factor  during  this  operation,  even  though  there  was  one 
injury.  Aircraft  flight  track  was  from  3  10  to  320’  during  this  jump,  w'hich  could 
have  exposed  the  jumpers  to  strong  cross  winds  if  the  winds  has  been  stronger. 
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Appendix  B 


List  of  Acronyms  ond  Symbols 


AAF 

- 

Army  Air  Field 

5WS  -  Fifth  Weather  Squadron  (AWS) 

AFGL 

AFSC 

— 

Air  Force  Geophysics 

Laboratory 

5WW  -  Fifth  Weather  Wing  (AWS) 

10th  Group  -  US  Army  10th  Special 

Air  Force  Systems  Command 

AGL 

- 

Above  Ground  Level 

ASL 

- 

Atmospheric  Sciences 

Laboratory  (Army) 

AWS 

- 

Air  Weather  Service 

Det 

- 

Detachment 

DZ 

GMT 

- 

Drop  Zone 

Cloud  Types 

- 

Greenwich  Mean  Time 

MEW 

- 

Mean  Effective  Wind 

Ci  -  Cirrus 

MSL 

- 

Mean  Sea  Level 

Cs  -  Cirrostratus 

NBC 

- 

Nuclear,  Biological,  Chemical 

Ac  -  Altocumulus 

NCAR 

- 

National  Center  for 

As  -  Altostratus 

Atmospheric  Research 

Sc  -  Stratocumulus 

NRC 

- 

Nuclear  Regulatory 

Commission  (US) 

Ns  -  Nimbostratus 

St  -  Stratus 

NWP 

* 

Numerical  Weather  Prediction 

PBL 

- 

Planetary  Boundary  Layer 

PE 

- 

Primitive  Equation(s) 

RAM 

- 

Random  Access  Memory 

SOWT 

- 

Special  Operations  Weather 

Team  (Air  Force) 

SP  -  Stability  Parameter 

TAP  -  Terminal  Aerodrome  Forecast 

TFRN  -  Terminal  Forecast  Reference 

Notebook 

TOT  -  Time  on  Target 

USAFETAC  -  IJS  Air  Force  Environmental  Technical  Applications  Center 
Z  -  "Zulu''  Time  (same  as  GMT) 


